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ABSTRACT 

We report on an analysis of the Extended Chandra Deep Field South (E-CDFS) region using archival 
data from the Very Large Array, with the goal of studying radio variability and transients at the sub- 
mJy level. The 49 epochs of E-CDFS observations at 1.4 GHz sample timescales from one day to 3 
months. We find that only a fraction (1%) of unresolved radio sources above 40 /iJy are variable at the 
4(j level. There is no evidence that the fractional variability changes along with the known transition 
of radio source populations below one milliJansky. Optical identifications of the sources show that 
the variable radio emission is associated with the central regions of an active galactic nucleus or a 
star-forming galaxy. After a detailed comparison of the efficacy of various source-finding algorithms, 
we use the best to carry out a transient search. No transients were found. This implies that the 
areal density of transients with peak flux density greater than 0.21 mJy is less than 0.37 deg~ 2 (at 
a confidence level of 95%). This result is approximately an order of magnitude below the transient 
rate measured at 5 GHz by Bower et al. (2007) but it is consistent with more recent upper limits 
from Frail et al. (2012). Our findings suggest that the radio sky at 1.4 GHz is relatively quiet. For 
multi-wavelength transient searches, such as the electromagnetic counterparts to gravitational waves, 
this frequency may be optimal for reducing the high background of false positives. 
Subject headings: catalogs — radio continuum: galaxies — surveys 



1. INTRODUCTION 

For more than four decades the largest science yield 
of variable and transient radio emission has come from 
single-dish radio telescopes, which have surveyed the sky 
for pulsed and periodic emission from pulsars and related 
phenomena on typical timescales of seconds to millisec- 
onds. In contrast, interferometric imaging surveys, which 
are best-suited to probe timescales ranging from seconds, 
hours and days, have only just begun. There are a limited 
number of surveys of these so-called long duration tran- 
sients and variables at frequencies of 1-10 GHz, e ach with 



al. 


2011 


Bell et al. 


2011 



Fully exploring this phase 
space is one of the main science drivers for a new genera- 
tion of synoptic radio imaging facilities, suc h as ASKAP 
dJohnston et al.|2008[ ) and Apertif/WSRT ( |Oosterloo et 



aT 2010[), that a re being built in the coming years (see 



aIso |Lazio et~aTp 009 ) 



Our knowledge of the variable GHz sky is especially 
lacking at sub-milliJansky flux density levels. For per- 
sistent sources there is a well-known flattening of the 
Euclidean-normalized radio source counts below about 1 
mJy, corresponding to a change in the radio source pop- 
ulations. This flattening is likely due to the emerging 
importance of star forming galaxies and low-luminosity 



( Condon 20071 [Seymour et al.| 


206^pmolcic et al. 2008 


Padovani 2011 


(Jondon et al. 


2012). Thus, while the 
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variability studies above milliJansky levels are dominated 
by rad io-loud AGNs with compact, flat-spect rum compo- 
nents (Sadler et al. 2006 Lovell et al. 20081, variability 
at these deeper flux density levels may probe new source 
populations. 

There are indications that the mJy transient sky is 
exciting. Over the last decade astronomers have de- 
tected transient decimetric emission from a variety of 
sources: transien t, bursting and pulsed radio emission 
from magnetars (j C ameron et al. 



2005 Camilo et aT'2006), short 



2005 Gaensler et al. 



ived radio afterglows 



ot short-duration gamma-ray bursts ( Fox et aLl|2005|) 



emission fro m (a transient) jet in a dwart nova ( Kording 



et al. 2008), a new population of sub-relativistic super 



novae ( Soderberg et al. 2010[), a mysterious popul ation 



of bursting radio sources H Hyman et al.|20 05, 2009 1, and 



erer et al. 


2011 


Cenko et al. 


2012) 



However, wi 



h few 



exceptions, most of what we know about the transient ra- 
dio sky has come via radio follow-up of objects identified 
by synoptic telescopes at optical, X-ray or gamma-ray 
wavelengths. Clearly, this titillating trove of (serendip- 
itous) discoveries call for systematic exploration of the 
decimetric sky on timescales of minutes and longer. 

A ready source of archival data for searching for sub- 
milliJansky transients and variables comes from deep 
continuum imaging surveys of the GHz radio sky un- 
dertaken to study the extragalactic radio source popula- 
tions. Many such surveys have been carried out, reaching 
noise levels of 4-10 /zJy and with arcsecond re solution 
(e.g. iSchinnerer et al.|[2007] |de Zotti et al.||2010| . In or- 
der to reach these deep flux density limits, it is standard 
to observe these fields for many epochs with a cadence 
that samples timescales of days, weeks and months. An 
added benefit is that these deep surveys are accompa- 



2 



nied by rich multi-wavelength continuum datasets and 
optical/infrared spectroscopic measurements. Thus the 
counterpart of any unusual variable or transient source 
can be readily identified and its redshift determined. 

In this paper we present a search for transients and 
variables at sub-milliJansky flux density levels using data 
taken as part of a deep radio continuum survey toward 
a region known as the Extended C handra Deep Field- 
South (E-CDFS; |MinCT^d][2008]). In p] we describe 
the original survey and our re-reduction of the data. In 
^3] we describe how the variability light curves for 599 
point sources were extracted. The transient search is 
described in S|4j The interpretation and implication of 
these results for radio source variability and transients is 
discussed in Sj5] 

2. OBSERVATIONS AND DATA REDUCTION 

The E-CDFS is an intensely studied region having a 
plethora of multi-wavelength data available (viz. X-ray, 
ultraviolet, optical, infrared and radio; see |Miller et al. " 



-27 : 20 : 00 



2008 



and references therei n) . Here, we use dat a from the 
radio survey undertaken by Miller et al. (2008j^J The ob- 
servations were made at a frequency of 1.4 GHz, using the 
National Radio Astronomy Observatory (NRAO) VLA 
in its A configuration (Project code AM 889). A hexag- 
onal grid of six pointings was made, with each pointing 
separated 12' from its nearest neighbor (see Table 111 and 
Figure [T]). All the observations were carried out in 2007 
between June 15 and September 23. Only a single point- 
ing was observed for each epoch. Each epoch was a 5-hr 
track centered on 03:30 LSTS There were a total of 49 
epochs, with a combined allocation of 245 hrs (see Ta- 
ble^. The total number of epochs for each pointing N ep 
is given m Table [TJ For more details abou t the s pecifics 



(2008) 



of the observational setup see Miller et al. 

The original purpose of these data was to average to- 
gether all pointings and epochs in order to create a deep 
(oVms = 5-8 /iJy) continuum image of the E-CDFS. In 
order to explore variability and to search for transients, 
we needed to work with the single epoch, single pointing 
images instead. Despite this, many of the data reduc- 
tion st eps that we followed were similar to |Miller et al.| 
(2008). We summarize the process here, pointing out 
slight differences. All calibration and imaging was car- 
ried out in the Astronomical Image Processing System 
(AIPS) packaged 

The calibratea uv data for each epoch were imaged 
and deconvolved separately. In order to image the entire 
field-of-view at full resolution, we created 37 different 
facets in a Ays-eye pattern. Each of the 37 facets was 
an image of 1024 2 pixels with a 0.5" pixel size offset 
from one another. There were another 23-25 outlier fields 
128 2 pixels each made of cataloged bright radio sources 
outside the primary beam of each pointing but within a 
2° radius. 

The AIPS task IMAGR was used to deconvolve each im- 
age to the rms noise level. In order to have a consistent 

5 In the pa st, radio observations o f this field have also been 
carried out b y Kcllcrman ct al. (2008}, [Afonso e t al. (2006}, and 
|Norris et al.| | |2UUbp . 

J This is of great benefit to variability studies, which are other- 
wise plagued by changes in the observing setup. 

7 http://www.aips.nrao.edu/ 
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Fig. 1. — Hexagonal grid of six pointings (crosses), each point- 
ing separated 12' from its nearest neighbor. The 34' X 34' (black 
square - ) region repr esents the extent of the final image from the 
|MilIer et al.| ( |2008| ) data release. 15'-circles (grey) corresponding 
to tne 507o beam attenuation in different pointings are also shown. 
The 60' X 60' background image is from NVSS. 



TABLE 1 
List of survey pointings 



Pointing ID 


R.A. (J2000.0) 


Dec. (J2000.0) 




ECDFS1 


03 33 22.25 


-27 48 30.0 


7 


ECDFS2 


03 32 55.12 


-27 38 03.0 


9 


ECDFS3 


03 32 00.88 


-27 38 03.0 


8 


ECDFS4 


03 31 33.75 


-27 48 30.0 


8 


ECDFS5 


03 32 00.88 


-27 58 57.0 


9 


ECDFS6 


03 32 55.12 


-27 58 57.0 


8 



Note. 



7Y cp is the number of epochs per pointing. 



set of images, we applied a Gaussian taper to the (30% 
level) of 100 kA and 70 kA in the u and v directions to 
the visibility data, and we restored the final images to 
a synthesized beam of 2.8" x 1.6" (position angle~ 0°). 
After deconvolution, the 37 facets for each epoch were 
combined using the AIPS task FLATN to form a single 
5120 2 pixel image 42.7' across. A correction for the at- 
tenuation from the primary beam was not applied at this 
stage in order that the images used in analysis had uni- 
form noise statistics across the entire image. The rms 
noise oy ms for each epoch is given in Table [2] 

3. VARIABILITY ANALYSIS 
The source catalog we used to investigate variability 



al.||2013 


) of 


Miller et al. ( 


2008) 



make a single dee p 34' x 34' image w ith a typical sensi- 
tivity of 7.4 /iJy. |Miller et al. (2013) identified sources 



using the AIPS task SAD down to 4cr, and then inspected 
the residual map to identify missed sources as well as ac- 
cepted sources which were poorly fit by SAD. These miss- 
ing sources were then added to the preliminary source 
list. Further flagging and follow-up was done in order to 
produce a modified source list in which all sources with 
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TABLE 2 
Observing Epochs 



Epoch 


Date 


Pointing oy ms 




UT 


fj, Jy bm^ 1 



1 


2007 


Jun. 


15 


H/CJJr bz 


26 


3 


2 


2007 


Jun. 


24 


hji^Db bo 


28 





3 


2007 


Jun. 


25 


H/L4Jr b4 


28 


6 


4 


2007 


Jul. 


01 


H/L4Jr bo 


25 


9 


5 


2007 


Jul. 


06 


ll/L4J^ bo 


29 


2 


6 


2007 


Till 

Jul. 


12 


H/L4Jr bl 


26 


6 


/ 


2007 


Jul. 


13 


TTr^TlTT'QO 

H/CJJr bz 


26 


6 


8 


2007 


Jul. 


14 


H/L4Jr bo 


25 


5 


y 


zuu * 


Jul. 


ID 


hiK^/Ur b4 


OR 
Zl) 


i 

i 


10 


2007 


Jul. 


16 


H/CJJr bo 


26 


8 


"1 1 

11 


2007 


Jul. 


17 


H/L4Jr bu 


26 


1 


12 


2007 


Till 

Jul. 


19 


H/CJJr bl 


26 





13 


2007 


Jul. 


20 


H/ODr bz 


34 


2 


14 


zuu * 


Jul. 


zl 


HjL4Jr bo 


OR 
ZD 


U 


"1 c 

15 


2007 


Jul. 




H/L4Jr b4 


25 


3 


16 


2007 


Jul. 


23 


H/L4Jr bo 


27 


8 


"1 T 


zuu * 


Jul. 


O A 
ZL 


hiK^/Ur bu 


0*7 
Z f 


O 


18 


2007 


Till 

Jul. 


26 


H/CJJr bo 


31 


8 


19 


2007 


Jul. 


27 


H/L4Jr bz 


27 


2 


o n 
zU 


zuu * 


Jul. 


OQ 
Zcy 


HjL4Jr bo 


n 
oU 


O 


zl 


2007 


Till 

Jul. 


29 


H/CJJr b4 


27 


1 


22 


2007 


Jul. 


30 


rL/L/JJr bl 


28 





O Q 
zo 


zuu * 


Aug 


no 
Uz 


IvbUr bO 


Ov 
Z f 


U 


o 1 


2007 


Aug 


03 


H/UJJr bl 


27 


4 


25 


2007 


Aug 


04 


17 r^Tl TP CO 

H/L4Jr bz 


26 


7 




zuu * 


Aug 


nt; 
UO 


rjbUr bo 


o X 

zo 


9 


Z 1 


onnv 
ZUU / 


Aug 


UO 


sliKjUr b4 


on 
zy 


1 


28 


2007 


Aug 


09 


T7r^T~lT7 1 Q C 

H/L4Jr bo 


31 


2 


29 


2007 


Aug 


10 


T7/^1T^iT7Q C 

H/ODr bu 


30 





■>n 
oil 


onnv 
ZUU / 


Aug 


1 1 


HjOUr bl 


■11 


n 
V) 


ol 


2007 


Aug 


13 


T7r^T~lT7 1 QO 

H/L4Jr bz 


31 


3 


32 


2007 


Aug 


14 


H/ODr bo 


3 1 


5 


QQ 

DO 


onnv 
ZUU / 


Aug 


id 


HjL/lJr b4 


■3 n 
oU 


Q 

o 


o4 


onn'7 
ZUU * 


Aug 




HjL4Jr bO 


4o 


•7 
f 


35 


2007 


Aug 


18 


ECDFS6 


28 


6 


36 


2007 


Aug 


21 


ECDFS1 


29 


1 


37 


2007 


Aug 


23 


ECDFS2 


30 


9 


38 


2007 


Aug 


25 


ECDFS3 


30 


8 


39 


2007 


Aug 


26 


ECDFS4 


29 


8 


40 


2007 


Aug 


28 


ECDFS5 


33 


9 


41 


2007 


Aug 


31 


ECDFS6 


29 


3 


42 


2007 


Sept 


. 06 


ECDFS2 


30 


9 


43 


2007 


Sept 


07 


ECDFS3 


29 


3 


44 


2007 


Sept 


. 08 


ECDFS4 


29 


3 


45 


2007 


Sept 


. 09 


ECDFS5 


28 


9 


46 


2007 


Sept 


10 


ECDFS6 


30 


2 


47 


2007 


Sept 


11 


ECDFS1 


30 


4 


48 


2007 


Sept 


12 


ECDFS5 


37 


9 


49 


2007 


Sept 


23 


ECDFS2 


31 


6 



Note. — List of the 49 epochs. Each epoch con- 
sisted of a 5-hr track centered on 03:30 LST. 



peak flux density greater than 5 times the local rms noise 
(i.e. 5a), were fit using the AIPS task JMFIT. Also, the 
effect of bandwidth smearing on sources within the six 
individual pointings was assessed using JMFIT, and the 
resolution information was thus preserved in the output 
catalog. Lastly, the sources in the DR2 catalog were 



(2008) 



compared with the catalog of Kellerman et al. 

The DR2 catalog contains almost twice as many 
sour ces (883 vs. 464) compar ed to the first data release 



(see Miller et al. 
data reduction. 



2013 
OTThe 



2008) owing to a more careful 



sources in the DR2 cata- 
log we created a point-source-only catalog of 736 objects 
used for exploring variability. With this careful approach 
outlined above, we expect the |Miller et al.| DR2 c atalo g 
to contain all real sources above 5a (however, see £4.1 1 



gating the variability of the sub-milliJansky population. 
In Sj4]we use both the DR2 image and its source catalog 
as a testbed for different source-finding algorithms. The 
signal-to-noise ratio (SNR) of sources in the DR2 catalog 
is shown in Figure [2j 



— PointSources 

— Extended Sources 

— All Sources 




Fig. 2. — A histogram of the signal-to-noise ratio (SNR) of the 
sources in the DR2 catalog. 

It is important to quantify the uncertainty in the peak 
flux density when analyzi ng variability of sources in dif- 
ferent epochs. Following Thyagarajan et al. (20111, we 
attribute this uncertainty primarily to a combination of 
foui^j causes: (i) local image rms noise, (ii) uncertainty 
in the primary beam, (iii) flux density calibration ampli- 
tude, (iv) pointing errors. Let us denote the measured 
peak flux density of a source in epoch i by such that 
the primary beam correction is b — b(9), where 9 is the 
angular distance of the source from the phase center. We 
wish to calculate the uncertainty in the quantity (fi/b). 
The image local noise (A/ 4 ) scales as A fi/b. If we denote 
the fractional uncertainty in the beam as eb and that in 
the flux density calibration amplitude by e c , then the 
corresponding errors scales as (fi/b)e. We adopt a value 
of 4% fo r e e , intermediate b etween the conservative esti- 
mate of Thyagarajan et al. 5%, and the one quoted by 



Ofek et al.| ( |201ip , 3%. "The typical pointing error (A0) 
of a VLA antenna is between 10"-20". The resultant 
uncertainty scales as (fi/b 2 )(—db/d9)A8. All these four 
error terms, added in quadrature, would give the total 
uncertainty, o^, in the peak flux density corrected for the 
primary beam attenuation. However, as shown below, 
the pointing-related error term is much smaller than the 
rest, and hence can be neglected. Thus, 



(1) 



We thus justify our use of the DR2 catalog for invest! 



Polynomial coefficients (and the associated error) that 

8 Far out in the primary beam, orthogonally polarized beams 
can be far offset on the sky, making amplitude calibration diffi- 
cult. For short observations, the effects of this "beam squint" can 
be much larger than all other uncertainties combined. However, 
these offsets tend to average out when observations are made over 
sufficiently long times. S ince the VLA's beam sq uint is oriented al- 
most exactly East- West | |Cotton"& : Pcrlcy 2010), our observational 
or averaging out tne ettects of 1 



setup is optimal for averaging i 



beam squint. 



4 



express the average angular dependence of b can be found 
in the AIPS task PBCDR, while measurements of the VLA 
be am power response to be yond the first null are given 



Cotton fc Perley | fl2010|). U sing the beam response 
profile trom |(Jotton fc Perley [ we can estimate the er- 
ror terms (i)-(iv) above for a typical source in the DR2 
catalog, having measured flux density of 300 ± 30 /i Jy. If 
the source lies at the half-power radius (9 = 15'), then 
these correction factors are about 10%, 4%, 4% and 1% 
of the primary-beam-corrected flux density (570 /iJy) re- 
spectively. The pointing-related uncertainty is thus neg- 
ligible. 

In light of Ofek et al. (|2011[), w e use two measures of 
variability (see also Scheers~|2011 ), the modulation index 



defined as the standard deviation divided by the mean 



11* 

7^ n^i^^"*^ 2 ' 



and the x 2 , 



N 



v2 _ y (ji - H 



(2) 



(3) 



where N is the number of epochs, and / is the mean flux 
density of the source over all the epochs considered for 
the variability analysis. In these two equations, the pri- 
mary beam correction is implicit in ji and /. x 2 gives a 
measure of the deviation from stochastic epoch-to-epoch 
fluctuations in the peak flux density, and we define "sig- 
nificant variability" beyond a leveHof 4a. The modula- 
tion index indicates the strength of variability, i.e. the 
fractional variation of the peak flux density. 

3.1. Single-Pointing Variability 

A variability analysis was carried out on each pointing 
in Table [l] separately. Peak flux densities were measured 
for all point sources in the DR2 catalog brighter than 40 
/zJy and within a 15' radius of the pointing centers (i.e. 
the 50% response radius of the primary beam of the VLA 
antennas; see also Figure]!]). This approach has the merit 
of being simple and robust. Since the angular distance of 
a source from its pointing center is constant, the accuracy 
of the correction for the primary beam attenuation b(9) 
is unimportant. The modulation index and \ 2 measures 
of variability are insensitive to a constant b(9). 

The limitation of this approach is that the resulting 
light curves are constructed for only 7-9 epochs. Many 
of our investigated point sources are found in multiple 
pointings and therefore light curves can be constructed 
with many more epochs, resulting in higher cadence over 
the full 100 days of observing . A full variability analysis 



of this kind is carried out in £ 3.2 



In Figure [3] we show variability plots (i.e 



X 2 vs. 



m) 

for the DR2~sources in each of the six pointings. There 
are approximately 175 sources per pointing. Some of 
the bright (> 3 mJy) sources show significant variabil- 
ity but with only low modulation indices (~10%). We 

9 For Gaussian noise, 4<r corresponds to a probability of about 
1/16,000, while the number of measurements in our variability 
analysis (several to tens of epochs multiplied by a few hundred 
sources) ranges from about 1,500 to 15,000. 



define strong variables as sources having m > 0.5 (i.e. 
higher than 50% fractional variability). Only two gen- 
uine variables were found in this single-pointing analysis; 
no strong variables were found (Table p] upper panel). 

3.2. Full Variability Analysis 

In order to undertake variability analysis using data 
from all the epochs, the single-epoch peak flux densities 
were required to be corrected for (i) bandwidth-smearing, 
and (ii) the primary beam response, as any given source 
will lie at different angular distance from the phase cen- 
ter in different pointings. We applied bandwidth smear- 
ing corr e ction from the approximation given in |Bridle fc] 
Schwabl ( p»89"l ) (Eqn. 13-19 therein). 

There were two choices for the primary beam profile 
- one derived as the empirical b eam profile for VLA- 
FIRST flThyagarajan et al.|[20TT| ), and the other found 
in the AlPS task PBC0R. To test which of these profiles 
best represented our data, we adopted the following ap- 
proach. We first normalized the peak flux densities of 
> 300 /xJy sources from all epochs using their peak flux 
density from the DR2 catalog and plotted them as a func- 
tion of distance from the pointing center. The resultant 
beam profile matched with the VLA-FIRST profile bet- 
ter than the one from PBCDR (to within 1%, but only for 
9 < 12'; scatter of 6.5%). Hence we used the former beam 
profile for our primary beam correcti on, b(9)\ the associ 



ated e rror (e&), was also taken from Thyagarajan et al. 
(2011 1. Thus, for a reliable all-epoch variability analysis. 



we restricted our search to the point sources in the DR2 
catalog which were located within 12' from the pointing 
centers of their respective epochs. This also appears to 
be the radius beyond which our bandwidth- smearing ap- 
proximation starts to break down. Thus, for example, a 
source located at a = 03 h 33 m 00 s and 5 = ~28°00'00" 
would be present in pointings 6, 5, and 1, but not in 2,3, 
and 4 because the separation between the source and the 
centers of pointings 2,3, and 4 is larger than 12 arcmin. 
F urth er, as we did with the single-pointing variability 
(§3TJ), we restricted our analysis to sources whose mean 
flux density was brighter than 40 /iJy. This full variabil- 
ity analysis was carried out on 599 point sources. 

Depending on the number of pointings in which a 
source is present, this analysis allowed us to exploit the 
higher cadence over the entire duration of the observ- 
ing program. The resulting light curves are now more 
densely sampled with 15-26 epochs, rather than the 7-9 
epochs for the single-pointing variability case. In Fig- 
ure^ we show variability plots for the DR2 sources, tak- 
ing into account all the epochs. 

Seven significant variables were found via this proce- 
dure, but no strong variables (m > 0.5). Both of the 
variables identified in the single-pointing analysis in §3.1| 
are also seen here. The results of the variability study 
are given in Table [3] and the light-curves for the signif- 
icant variables are shown in Fi gure [5| We can com pare 
our variability criteria with the CarillT et al.| ( |2003[) mea- 
sure for variability, i.e. V c = (Si — S2)/S, where d'i and 
5*2 are the maximum and minimum flux densities respec- 
tively observed, and S is their mean. This measure of 
variability for the seven variable sources found in this 
work are listed in Table |3] 

In order to undertake multiwavelength identifications 
we had to align the reference frames of all the data sets. 
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Fig. 3. — Variability plot for the point sources in the Miller ct al. DR2 catalog, shown separately for all pointings. The peak flux density 
is denoted by the symbol size. The red dashed line represents the 4cr level for the appropriate number of degrees of freedom (one less the 
number of epochs in each pointing) for each pointing. The number of epochs in each pointing lies between 7 and 9. 



TABLE 3 

Variables among Miller et al. DR2 sources 



ID 


Q J2000 


<5,72000 


/ 


m 


V c 




z log L R 






Mr 


Energy 




(h, m, s) 


(V,") 


0*Jy) 




(%) 






(cgs) 








(mag) 


Source 


Single Pointing Analysis 


337 


03 32 18.03 


-27 47 18.8 


558 ± 14 


0.26 


76 


8 
















621 


03 33 15.00 


-27 51 51.3 


663 ± 22 


0.23 


93 


9 
















Full Variability Analysis 


67 


03 31 27.07 


-27 44 09.9 


136 ± 9 


0.50 


181 


16 


1.005 


30.6 


+0.23 


> 


-3.8 


-21.04 


SF+AGN 


190 


03 31 52.13 


-27 39 26.6 


891 ± 12 


0.15 


48 


16 


2.296 


32.1 


-0.35 


> 


-1.3 


-22.07 


AGN 


239 


03 32 00.85 


-27 35 57.1 


1822 ± 27 


0.12 


38 


8 


0.266 


30.5 


-0.30 


> 


-2.3 


-17.26 


SF+AGN 


297 


03 32 11.66 


-27 37 26.3 


3600 ± 36 


0.11 


41 


17 


0.605 


31.5 


+0.89 




-0.9 


-23.98 


AGN 


337 


03 32 18.03 


-27 47 18.8 


494 ± 11 


0.26 


95 


16 


0.734 


30.8 


-0.08 


> 


-3.4 


-22.69 


SF+AGN 


621 


03 33 15.00 


-27 51 51.3 


497 ± 10 


0.20 


70 


15 


1.107 


31.2 


+0.24 


> 


-3.4 


-20.71 


SF+AGN 


628 


03 33 16.74 


-27 56 30.4 


1341 ± 16 


0.11 


38 


15 


0.685 


31.2 


-0.40 




-2.2 


-21.21 


SF+AGN 



Notes 

the primary beam and ba ndwidth smear i ng 
Redshift z is according to |Bonzini et al.l J 201J2 



(3) 



Vp is the Carilli et al. l 



Komatsu ct 



. (2) / refers to the mean flux density corr ected for 
( |2003| > variability criteria as described in § |3.2| (4) 
or |Treister et al.| ^UTOf . The redshifts of ID 67, ID 239 and ID 297 as per 

the COMBO-17 survey catalog fWolt et al 

(erg s — 1 Hz -1 ) is Lr = 47rd^//(l + z), where the luminosity distance d[ assumes cosmological parameters from 
[aL] l |2011| . (5) The spectral indices (S oc v a ) between 1.4 GHz and 5.5 GHz and between 24 fim and 70 fim are tabulated as 
ctpi and qir respectively. In the absence of a 70 /im counterpart, a 3cr upper limit to the flux density at this wavelength is 
considered. Note that the 1.4 GHz and 5.5 GHz measurements are non-simultaneous. (6) The absolute R-band magnitude Mr 
has been calculated using the redshift z and the apparent magnitude from the NASA/IPAC Extragalactic Database. (7) The 
energy sources within the galaxy as indicated by t he ra dio and mid-to-far-infrared properties are listed in the last column. SF: 
star-formation, AGN: active galactic nucleus (see §331. 



G 



2.5 
2.0 
~1.5h 
"x 1.0 
0.5h 



• 

• 

• 




• • 0.1 mjy 
° « 0.3 mjy 
o o l mjy 
o o 3 mjy 
O O 10 mjy 


• 




• 




o 

8 \o 

CD ..Pfl-. „ 







0.2 0.4 0.6 O.i 
m 



1.0 



1.2 



1.4 



Fig. 4. — x normalized by its value at the 4c l evel, plotted 
against the modulation index, m, for sources in the |Miller et al.| 
DR2 catalog using peak fluxes from all pointings, and corrected 
with empirically-derived beam attenuation profile. The 4(7 level is 
different for different sources, depending on the number of epochs in 
which they are present. The mean peak flux density is denoted by 
the symbol size. Filled circles indicate significant variables (lying 
above a normalized \ 2 of unity). 



The radio and optical source positions were brought to 
the same reference frame by calculating the radio-source 
position offsets with respect to Hubble Space Telescope 
source positions from GOODS-S ( Giavalisco et al.||2004 | 
and the GEMS (Rix et al. 2004). Optical counterparts 
(from these two HST catalogs) were searched toward ra- 
dio sources within 1". Only point-like sources having 
a single counterpart were chosen. A histogram of the 
offsets of these counterparts in right ascension (Aa = 
aradio - ^optical) and declination (AS = S iad i - Optical) 
was then computed to find the most-likely offset (peak of 
the histogram) . The associated error was taken to be the 
standard deviation about this most-likely offset added in 
quadrature with the cataloged mea n positional error of 



the radio source. In addition to the Miller et al. (20131 
DR2 catalog, we repeated this procedure fo r other radio- 
sou rce ca t alogs in the E - CDFS regi on, viz. |Kellerman et 



al. (2008), Norris et al. (2006), and Afonso et al. (2006) 



using ^5<7 sources. The computed radio vs. optical posi 
tional offsets along with the mean positional errors listed 
in the respective radio catalogs are shown in Table [4] 



TABLE 4 

Radio Positional Offsets wrt. Optical HST Catalogs 







(1) 


(2) 


(3) 


(4) 








~ 0.1 


0.37 


0.31 


b 






"S 


~ 0.1 


0.57 


0.58 


b 




GEMS 


Aa 


0.18±0.31 


0.15 ± 0.45 


-0.09±0.41 


-0 


17 ± 0.40 a 


AS 


-0.32±0.32 


-0.34 ± 0.53 


-0.20±0.44 





11 ± 0.30 a 


GOODS 


Aa 


-0.20±0.23 


0.15 ± 0.42 a 


0.03 ± 0.44 a 


-0 


10 ± 0.39 a 


AJ 


0.22±0.28 


-0.18 ± 0.37 a 


0.03 ± 0.44 a 


-0 


18 ± 0.33 a 
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Fig. 5. — Light curves of the variable sources from [Miller et ah] 
DR2 catalog listed in Table [3] Filled symbols represent the flux 
densities considered for variability analysis (i.e. where 8 < 12'). 
Flux densities in epochs where a source lies within the 50% power 
circle of the beam (i.e. where 9 < 15') are plotted for reference as 
open symbols. The error bars take into account the background 
rms, primary beam correction, and bandwidth smearing (no taper). 

3.3. Notes on Variables Found 

An important question that we can address is whether 
the variability at sub-milliJansky levels is dominated 
by normal star-forming galaxies or by AGN. For AGN- 
dominated samples above 1 mjy, variability at frequen- 
cies of a few GHz or below is thought to be dominated by 
propagation effects (i.e. refractive interstellar scintilla- 



tion) and not by intrinsic changes in t he source ( Gaensler 
fc Hunstead|2000| |Ofek fc Frail|201l| ). For a disk galaxy 
we expect there to be steady emission from diffuse syn- 
chrotron emission and the sum of all supernovae. Intrin- 
sic variability can be induced by (i) a stellar explosion 
(supernova, low-luminosity gamma-ray burst), and (ii) 
nuclear radio emission (AGN). Mapping the radio emis- 
sion to the center of the galaxy would favor (ii) whereas 
if the radio emission is mapped to the disk then (i) is 
favored. 

To this end, we overplotted radio source positions 
of the seven variables given in Table [3] on HST image 
cutouts from GOODS-S and the GEMS projects (Fig- 
ure [6]) . All of the radio sources have an optical coun- 
terpart on these HST images. The redshifts of these 
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Fig. 6.— 2.2" x 2.5" GEMS (all sources except ID 337) and GOODS-S (ID 33 7) F606W ACS-WF C image cutouts of variable sources 
found in this work. The red error ellipses denote the radio source positions from [Miller et aT| ( |2013} , shift-corrected to the HST source 
positions. All the positional uncertainties are lcr (see § |3.2| and Table H|. 



objects vary from 0.3 to 2.3. Thus, the post-offset radio 
source-position uncertainties lie between 1.4 and 2.7 kpc. 
Within most error ellipses, there is a centrally compact 
source, suggesting that this is the source of the variable 
emission. Light-curves of the variable sources are given 
in Figure [5] Detailed notes on each of these objects are 
given below, and key physical parameters are listed in Ta- 
ble |U All of the variable sources have luminosities in ex- 
cess of 10 30 erg s _1 Hz -1 , where the luminosity functions 



of AG Ns and star- forming galaxies intersect (Condon 
2002 1 . Most of the variable sources have a roughly flat 



spectrum between 1.4 GHz and 5.5 GHz suggesting the 
presence of AGN. Note that the spectral indices between 
1.4 GHz and 5.5 GHz are based on non-simultaneous 
measurements having different resolutions. Addition- 
ally, mid-infrared colors, far-infrared spectral indices, 
and mid-to-far-infrared luminosities indicate that most 
of the variable sources reside in star-forming galaxies. 
In Figure [7] we show the mid-infrared color-color dia- 
gram for the variable sources using Spitze r/IRAC pho 



tome try fro m the SIMPLE survey catalog (Damen et al. 
2011| . The |Donley et aT] ( |2012| AGN selection region 



is overplottecT lor comparison, th e colors of the radio 
sources from the AEGIS20 sample ( |Willner et al]|2012 1 
are shown along with those of three template spectral en- 
ergy distributions (SEDs) — an elliptical galaxy, an Sbc 
galaxy, and an AGN — from |Assef et al.| ( |2oTo| ). This 
figure along with the far-infrared flux densities identify 
two variables with AGNs and the rest as star-forming 
galaxies. 

Summarizing, the high-resolution optical images to- 
gether with photometric information from radio and mid- 
to-far-infrared suggest that variability arises from the 



central regions of an AGN or star-forming galaxy. 
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Fig. 7. — Spitzer/IRAC color-color diagram for the mid-infrared 
counterparts of the variable sources (fill ed circles). For com par- 
ison, the radio sources from AEGIS20 jWillner et al.||2012| are 
also shown (unfilled circles). The AGN selection region m the 
upper - right corner, bounded by the dashed lines, is from |DonIey| 
|et al.| 1^2012 ) . Curves represent the colors of three template spec- 
tral energy distributions (E: el liptical galaxy; Sbc: spiral galaxy; 
AGN) from |Assef et al] pOlO) as redshift increases from to 3. 
The redshift-dependent color-coding is — blue: z ^ 0.5, green: 
0.5 < z < 1.1 and red: z > 1.1. 



ID 67: This source has a mean 1.4 GHz flux density 
of 136 ± 9 /iJy and shows the strongest variation among 
all seven variables. Two flaring bursts are evident from 
the light curve (Figure 5]), one of them lasting for about 
25 days, and the flux density of the other increasing by 



8 



at least a factor of 2.5 in 12 days. This source also ap- 
pears in the Kellerman et al. 1.4 GHz catalog, where its 
flux dens ity is 90 ± 16 juJy. T he photometric redshift of 
z=1.005 " 



nosity of 



Bonzini et al.|2012 ) implies a mean radio lumi- 
30.6 erg s _1 



og Lr = 30.6 erg s Hz . The radio spectral 
index (defined as a where S„ oc u a ) be tween 1.4 GHz and 
5.5 GHz is +0.23 ( |Huynh et aL|2012|) The morphology 
of the host galaxy as seen from the HST image appears 
to be that of a bright, compact nucleus surrounded by 
faint extended structure. The half-light radius of the 
galaxy according to the GALFIT parameters from the 
GEMS catalog is 620 ± 10 pc. The K-correctecH mid- 
to-far-infrared spec tral luminosities derived fr om the FI- 
DEL and GOODS ( |Magnelli et al.|2009|[20lT| , and SIM- 
PLE Spitzer surveys are vL„(5.8 /jm) = 6.6 x 10 8 L Q 
and vL v (2A /im) < 5.4 x 10 10 L^ . By c omparing these 



quantities with the Chary & Elbaz (2001 1 template SEDs 
(see Figure 4 of that paper) and from the Spitzer/IRAC 
color-color diagram (Figure [7]), we interpret that the host 
galaxy is star-forming. 

Taken together, the radio luminosity, radio spectral 
index, optical morphology and mid-to-far-infrared flux 
densities argue that ID 67 is a star-forming galaxy har- 
boring a low-luminosity AGN. 

ID 190: This source has a mean 1.4 GHz flux density 
of 891 ± 12 /iJy. The light curve appears to fluctuate be- 
tween high and low flux density states on a timescale of 
order 20 d ays. This sour ce is also found in the Kellerman 
et al. and Norris et al.| catalogs where the flux density 
is 970 ± 30 and 810 ± 19 /xJy respectively, and a spec- 
tral index of -0.6 bet ween 1.4 and 4.8 GHz. The 5.5 
GHz flux density from |Huynh et aL is 555 ± 17 /zJy, for 
which we derive a spectral index of —0.35. T he photo- 
metric redshift of z=2.296 ( Bonzini et al.|20T2 ) implies a 
mean radio luminosity of log Lr = 32.1 (cgs). The HST 
image of this galaxy as seen in Figure [6] shows that the 
radio source is offset from the brightest emission in the 
field. On t he basis of its departure from the radio-FIR 
correlation, |Norris et al.| classify this source as an AGN. 
The radio and far-infrared spectral indices (Table |3| and 
mid-infrared colors (Figure [7]) are consistent with this 
identification. Note that between 24 /zm and 70 /im, 
a spectral i ndex greater t han —1.5 is representative of 
AGNs (e.g. lCondon"|[2002l ). 

ID 239: This source has a mean 1.4 GHz flux density 
of 1822 ± 2 7 /xJ y. This sourc e is also found in the |Keller^1 
man et al. and |Norris et al. catalogs where the flux den- 
sity is 2030 ±43 and 1640 ±20 £tJy respectively. Its light- 
curve shows a gradual increase in flux density, over a pe- 
riod of about 50 days, followed by an equally gradual de- 
cline. The radio spectral index between 1.4 GHz and 5.5 
GHz is —0.30. This slope is more consistent with AGNs 
than star-forming galaxies whic h are expected to have 
a ~ -0.8. The redshift z=0.266 ( |Treister et al.|2009[ ) im- 
plies a mean radio luminosity of log Lr = 30.5 (cgs). The 
HST cutout (Figure [6]) reveals just a bright nucleus sur- 
rounded by faint extended structure, the half-light radius 
being 3.2 ± 0.5 kpc. The K-corrected mid-to-far-infrared 
spectral luminosities, vL v (5.8 /im) = 3.6 x 10 s L Q and 
vL„(2A fim) < 1.0 x 10 9 L , and the mid-infrared colors 

10 For the 5.8 (im and 24 /im luminosities, the K-correction has 
been applied based on the 8 fim and 70 /im flux densities. In the 
absence of a 70 (im detection, the 3<r upper limit has been used. 



(Figure [7]) show that the host is a star- forming galaxy. 

Taken together, the radio luminosity, radio spectral 
index, optical morphology and mid-to-far-infrared flux 
densities argue that ID 239 is a star-forming galaxy har- 
boring a low-luminosity AGN. 

ID 297: This source has a mean 1 . 4 GH z flux den- 
sity of 3.60 ± 0.04 mJy. Huynh et al. (20121 measure a 
5.5 GHz flux density of 12.25 mJy, implying a spectral 
index of + 0.89. One the basis o f this steep positive spec- 
tral index, Huynh et al. (2012) suggest that this is part 
of a class of Gigahertz Peaked Spectrum (GPS) s ources, 



thought to be a young AGN. The redshift z=0.605 ( Treis 



ter et al. [2009 ) implies a mean radio luminosity of log 
Lr = 31.5 (cgs). ID 297 appears to be a stochastically 
varying source with no specific trend in its light-curve 
(Figure [5]). The radio position is consistent with a bright, 
unresolved HST source. On the basis of departure from 
radio-FIR correlation |Norris et al.| classify this source 
as an AGN. The mid-infrared colors (Figure [7]) and far- 
infrared spectral index (Table [3]) are consistent with this 
identification. 

ID 337: This source has a mean 1.4 GHz flux density 
of 494 ±11 /iJy. From its radio light-curve, ID 337 ap- 
pears to have repeated outbursts roughly every 25d. The 
flux density of this sou rce in the|Miller et al.{ Kellerman 
et al. Norris et al. and Afonso et al.| catalogs is 439 ± 8, 



524 ± 14, 380 ± 16, and 404 ±34 /iJy respectively, and a 
spectral index of —0.2 between 1.4 and 4.8 GHz. The 5.5 



GHz flux density from |Huynh et al. is 443 ± 20 \i Jy, for 
which we derive a spectral index of —0.08. Th e photo- 
metric redshift of z=0.734 ( Vanzella et al.|[2008 ) implies 
a mean radio luminosity of log Lr = 30.8 (cgs). 

The HST image (Figure[6]) shows that the optical coun- 
terpart to the radio source is actually associated with 
the fainter (or more extinctcd) galaxy among a group of 
two closely separated galaxies. This has lead to some 
confusion over the correct optical identification . VL A 



observation s car ried out in 19 99-2001 and 2007 (Keller- 

respectively) suggest that 



man et al. and Miller et al 



the radio counterpart is the fainter galaxy. ATCA ob- 
servations car ried out in 2003-2004 ( |Norris et al. and 
Afonso et al.| ) suggest the bright galaxy as the opti 



cal counterpart. Our optical-radio frame tie summa- 
rized in Table [4] supports the fainter optical source as 
the likely radio counterpart. The spectral index and 
radio luminosity arg ue that ID 337 is a low-luminosity 
AGN. |Afonso et al. state that the bright radio source 
is a luminous star-forming galaxy, possibly part of a 
merging system. Comparison of the K-corrected spec- 
tral luminosities, z/L„ (5.8 /mi) = 4.0 x 10 9 L and 
10 10 



vLjJ^A jim) < 1.1 



L^, with the 



Char> 
(2007) 



fc Elbaz 
together 



(2001| SED templates and Desai et _al i 
with the mid-infrared colors (Figure[7| advocate the star 
forming nature of the host galaxy. 

ID 621: This source has a mean 1.4 GHz flux density 
of 497 ± 10 /iJy. The first epoch in its light-curve reveals 
an initial brightening of the source followed by a decline 
and subsequent small-amplitude variations. The maxi- 
mum flux density is close to 4 
density of about 300 /xJy. The 



Miller et al. 


Kellerman et 


ist the flux c 


ensity of this 



aL} and |Norris et al.| catalogs '. 
source as 494±10, 565±17, 450±18 /iJy respectively, and 
a spectral index of —0.1 between 1.4 and 4.8 GHz. The 
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5.5 GHz flux density from |Huynh et aL is 689 ± 16 /iJy, 
for which we derive a spect ral index of 0.24. Th e photo- 



metric redshift of z=1.107 (Bonzini et al. 20121 implies 
a mean radio luminosity of log Lr = 31.2 (cgs). 

The HST image reveals a bright nucleus surrounded by 
diffuse emission. The half-light radius is ^3 kpc. The ra- 
dio position is consistent with the nuclear source. Taken 
together, the radio luminosity, spectral index and optical 
morphology argue that ID 621 is an AGN. Additionally, 
the K-corrected 5.8 /jm luminosity of 5.0 x 10 9 L , the 
24 /im luminosity upper limit of 6.8 x 10 10 L Q , and the 
mid-infrared colors (Figure [7| suggest that this is also a 
star-forming galaxy. 

ID 628: This source has a mean 1.4 GHz flux density 
of 1.34 ± 0.02 mJy. The light-curve indicates a steady 
increase of the quiescent emission to maxi mum flux den- 
sity, followed by a s teady decline . The |Miller et al.[ 
|Kellerman et al.[ and |Norris et al.| catalogs list the flux 
density of this source as 1.07 ± 0.02, 1.33 ± 0. 03, and 
0.90 ± 0.02 /iJy. The 5.5 GHz flux density from |Huynh| 



et al. is 0.78 ± 0.02 /xJy, for which we derive a spectral 
index of —0. 40. T he photometric redshift of z=0.685 



( Norris et al. 2006 ) implies a mean radio luminosity of 
log L R = 31.2 (cgs). 

The HST image reveals a bright nucleus surrounded by 
disk-like diffuse emission, the half-light radius being ^10 
kpc. The radio position is consistent with the nuclear 
sour ce. On the ba sis of departure from radio-FIR correla- 
tion [NorrisjeFar] classify this source as an AGN, which is 
consistent with our flat radio spectral index and radio lu- 
minosity. Additionally, the K-corrected 5.8 /mi luminos- 
ity of 3.0 xlO 9 Lq, the 24 \xm luminosity of 3.7xl0 10 L Q , 
and the mid-infrared colors (Figure [7]) suggest that this 
is a normal star-forming or a starburst galaxy. Taken to- 
gether, the radio and mid-to-far-infrared properties along 
with the optical morphology argue that ID 628 is an AGN 
embedded within a star-forming galaxy. 

4. TRANSIENT SEARCH 

For our transient search, we are interested in identi- 
fying those point sources which show up above the flux 
density limit for a short amount of time (correspond- 
ing to one or more epochs depending on the cadence), 
and remain below the limit in all other epochs. It is 
therefore important to reliably distinguish noise from 
real transients (which will determine how well we can 
reject false positives), and to avoid the rejection of tran- 
sients as noise (reducing the number of true negatives). 
Thus, characterizing the effectiveness of source-finding 
algorithms in terms of reliability and completeness is 
crucial. Another motivation for characterizing source- 
finding algorithms is to find the optimum parameter val- 
ues for use in transient searches. 

In the following subsections we begin by the te stin g 
the efficacy of existing source-finding algorithms (^4.1). 
We then apply some of th e be tter-performing algorithms 
on the E-CDFS epochs (JO). 



4.1. Efficacy of Source- finding Algorithms 

Recently, quantitative tests have been carried out on 
the reli ability and completeness of source-finding algo- 



and IMSAD in MIRIAD, Aegean, SExtractor, and Selavy, 
these studies found that sf ind and Aegean produce the 
most reliable cata logs 



T he an alysis of Huynh et al. (2011) and Hancock et 



al. d20l2) was carried out on two simulated data sets 
JT) The ASKAP simulation is a 4° x 4° image of a full 
continuum observation with critically-sampled beams in 
the 6-km ASKAP configuration. Its pixel-scale is 2.75" 
and the rms noise is about 35 mJy, which varies across 
the field. 16 idealized beams one degree apart, spaced 
in a rectangular grid, mimic the effect of the phased- 
array feed. The image contains ~7.7 million sources hav- 
ing flux densities greater than 1 /j jy from the S 3 -semi - 
empirical extragalactic simulation ( Wilman et al.pblO ). 
(ii) For the Hancock et al. simulation, a sky image was 
created as a 4801 2 pixel image 8° across with a 6" pix- 
els sampling a 30" beam and an rms noise of 25 /iJy. 
Sources were injected at random positions with angular 
sizes (with random position angles) from 0" to 52" and 
with source number counts distributed with peak flux 
densities as N(S)cx S~ 2 - 3 such that 15,000 sources having 
fluxes densities > la are present in the image. 

With our E-CDFS dataset we are able to carry out 
a similar analysis using real data with all its attendant 
residual calibration and imaging errors. A comparison of 
real and simulated data could be informative. While sim- 
ulations are useful in determining which source-finding 
algorithm work best in general, they do not explore the 
parameter space of the algorithm thoroughly. Hence they 
may not provide optimum parameter values for a tran- 
sient search on a specific dataset. The deep field of the 
E-CDFS is well suited for this comparison. The 4096 2 
pixel image is 34' across with a 0.5" pixels sampling a 
synthesized beam of 2.8" x 1.6" (position angle~ 0°) and 
an rms noise 7.4 ^Jy. As outlined in ^J3J great care was 
taken in constructing the DR2 source catalog so we can 
be assured of its completeness and reliability (see also 
Figure 2]). 

In what follows, we will use the DR2 catalog and the 
deep E-CDFS image to test various source-finding algo- 
rithms for completeness and reliability for different input 
parameters. We use the terms "real sources" and "false 
sources" as being those sources present in the DR2 cat- 
alog, and those that are not, respective!}^] Whether a 
source detected by a source-finding algorithm has a coun- 
terpart in the DR2 catalog is determined by searching for 
DR2 sources within 1" from the source position. Follow- 



ing Huynh et al. (2011 1, we define two additional terms, 
"completeness" , as the fraction of real sources detected 
by a source-finding algorithm, and "reliability", as the 
fraction of detected sources which are real. Note that 
in Figures [5J [9] and [T2j we plot the complet eness and 
reli ability with in contiguous SNR bin s, unlike [Huynh et 



rithms (Huynh et al. 2011 Hancock et al. 2012). Of 



these many publicly available software packages, sf ind 



al. ( 2011[ ) and |Hancock et aL (2012) where, for a given 
SINK, the plots represent the completeness and reliabil- 
ity for sources greater than or equal to that SNR. For 
source-finding algorithms employing a probabilistic ap- 
proach of drawing a pixel from the background and thus 

11 Even though the DR2 catalog was constructed with great 
care, it is likely that ~1 genuine source was missed and a handful 
of spurious sources added (perhaps not truly spurious, but SAD 
sources at 4<r bumped up to 5<r by the JMFIT task in AIPS). Recall 
that according to Eddington bias, more sub-5<r sources get bumped 
up than 5<r get bumped down. 
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calculating false-detection rate (FDR), FDR + Reliabil- 
ity = 100%. 

4.1.1. sfind (MIRIAD) 

In its default mode, sf inc^] incorporates a statisti- 
cally robust method for detecting source pixels, called 
"False Discovery Rate" , or FDR. In the FDR algorithm, 
detected sources are drawn from a distribution of pixels 
with a robustly known chance of being falsely drawn from 
the background. Contiguous, monotonically decreasing 
adjacent pixels from the FDR-selected ones, are used for 
fitting 2-D elliptical Gaussians to the sources. Thus, the 
fraction of expected false sources is more reliably deter- 
mined than in sigma-clipping criteria (see metho ds be- 



also f 4.1.3 for a note on the efficacy of SAD). These find- 



low). Details of the FDR method can be found in Hop 
kins et al. (2001). The run-time for sfind searching tor 
sources down to 5a in the DR2 image is about 20 sec- 
onds. For all the tests carried out on sfind, MIRIAD 
version 4.2.3 (optimized for CARMA; CVS Revision 1.11, 
2011/04/26) was used. 

Here, we explore the completeness and reliability of 
sfind by tweaking the two relevant parameters: (i) 
alpha (a), the percentage of probable background pixels 
that can be accepted in the analysis, and (ii) rmsbox, the 
size of the smoothing box used for estimating the back- 
ground and the standard deviation of the image. Table [5] 
lists the parameter values tested. 

In general, completeness of sfind increases, and its 
reliability decreases (Table ^ with increasing a, as ex- 
pected for an FDR algorithm. Also, both these quantities 
increase with rmsbox. The rise in completeness and reli- 
ability is precipitous (between 5% and 20%) as rmsbox is 
increased from 5 to 10 beamwidths, after which it flattens 
off. There is a slight decrease in the reliability as rmsbox 
is increased from 20 to 50 beamwidths. 10 beamwidths 
can then be interpreted as the minimum box size for 
determining the background rms noise without signifi- 
cant contribution from the sources themselves, whereas 
50 beamwidths would correspond to the size where the 
calculated rms starts deviating significantly from the true 
local rms. However, in the case of bright sources in the 
field, these rmsbox limits might be somewhat larger. 
Completeness and reliability are ~88% for a=10 and 
rmsbox— 20, but better reliability (by a few percent) can 
be obtained at an equal expense of completeness by using 
a=5 or 2. Further decrease in a substantially decreases 
the completeness without any significant improvement 
in the reliability. Thus, we determine the optimum val- 
ues for the input parameters to be 5 ^ a ^ 10 and 
10 ^ rmsbox ^ 50 beamwidths. For the FDR algorithm, 
we expect the reliability to be 1 — a. However, sfind 
reliability is lesser than this expected value by a few per- 
cent, due to the acceptance of sidelobes of bright sources 
and fitting of extended sources with multiple elongated 
and overlapping (unphysical) components. To some ex- 
tent, the acceptance of these false sources can be reduced 
by choosing a lower value of a. Some examples of the 
components fit by sf ind to different kinds of sources in 
different environments in the DR2 image are shown in 
Figure [6] In general, compared to other source-finding 
algorithms excepting SAD, sfind has a better reliability 
for a given completeness, and vice-versa (Tables [§|9| see 

12 www.atnf.csiro.au/computing/ software /miriad / doc /sfind. html 



ings are consistent with the tests carried out on simulated 
images ( Huynh et al.|[20TT ). In the ASKAP simulation, 
the reliability of the sf ind catalog is rather flat with re- 
spect to a, whereas, for the Han cock et al.| simulation it 
decreases by severa l percent as a increases from 0.1 to 10 
( Huynh et al.|20lT |. In the latter, rmsbox— 20 is found to 
give slightly better reliability than rmsbox=20. For the 
~3-10ct sources in these simulations, the completeness 
for a— 5 is greater than that for a=0.1 by 5-10%. They 
do not explore the completeness for a=10 and rmsbox 
other than 10 beamwidths. 

Figure J5] shows the completeness of the sfind cata- 
log as a function of the SNR of the detected sources for 
different values of the input parameter a. The reduc- 
tion in completeness beyond SNR~20 is due to missed or 
badly-fit components of extended sources. The complete- 
ness for optimum values of input parameters is shown in 
the upper panel of Figure |9j which shows that sfind 
reaches 100% completeness at a much lower SNR than 
other algorithms. The lower panel of Figure M shows the 
reliability. Curiously, there is a dip in the reliability be- 
tween SNRs of 6 and 10. Point sources strewn across 
the DR2 image but primarily located near its edges (in 
regions of increased rms; see the lower panel of Figure [10] 
for example), which are absent in the DR2 catalog and 
detected by sfind, are responsible for this reduced relia- 
bility. Some of these sources get rejected from the sfind 
catalog when a lower value of a or rmsbox is used, indi- 
cating that they are either not genuine sources on the sky 
or are sources at a lower SNR. Indeed, other algorithms 
find some of these sources to be at a much reduced SNR 
between 3 and 5. The results of (Hancock et al. 20121 
also show a dip in th e sf ind reliability, although centered 
on SNR^15. Figure |TT| compares the SNR of sources de- 
tected by sfind, with those of the counterparts from the 
DR2 catalog. The largely increased SNR reported for 
sources at or below 6a in the DR2 catalog, and largely 
reduced SNR for sources above ~100cr, is anomalous. Al- 
though for SNR above 100 the discrepancy is likely to be 
due to overlapping components fit to extended sources, 
all the above observations indicate that the SNR reported 
by sfind is somewhat different from that reported by 
other algorithms. Sources with 6<SNR<10 in the sfind 
catalog clearly have over-estimatated (upto 100%) peak 
flux densities and different dimensions than what would 
be expected by inspecting the cutouts of these sources. 
Thus, fitting of sources rather than differently-reported 
rms is the cause of the SNR descrepancy, at least at these 
low SNRs. The solution of this issue is possibly in cor- 
rectly setting the fdrpeak and psfsize input parameters 
which allow reasonable measurements of sources close to 
the threshold. Since there are not many sources in our 
data having SNR >70 (Figure [2]), the suggested trend 
in completeness, reliability and measured SNR in this 
domain should be treated with caution. 

4.1.2. S Extractor 

SExtracto^] is a source-finding program widely used 
in optical astronomy, and is particularly oriented towards 
the reduction of large-scale galaxy-survey data, as well 

13 www.astromatic.net / software / sextractor 
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as sparsely- to moderately-crowded stellar fields. It ana- 
lyzes the image in two passes such that, in the hrst pass, a 
background map is made, and in the second, background- 
subtraction, filtering and thresholding is done on-thc-fly. 
Detected sources are then deblended and CLEANed be- 
fore performing photometry. The run-time for SExtrac- 
tor searching for sources down to 5a in the DR2 image 
is about 2 seconds. For all the tests carried out here, the 
latest release of SExtractor (version 2.8.6, 09-Apr-2009) 
has been used. 

We adopted a strategy of searching for sources down to 
Nct (N=3,5,6,7,10, as set using the DETECT .THRESH 
parameter in SExtractor) and then selecting the 
ones ^5cr prior to the comparison with the DR2 
sources. Two input parameters (apart from possibly 
DEBLEND.NTHRESH, BACK.FILTTHRESH and 
C LEAN _P ARAM, which we have not tested) are ex- 
pected to affect the completeness and reliability of the 
SExtractor catalog. First, the mesh-size, which deter- 
mines the size of the box used for background rms es- 
timation, is specified by BACKSIZE, and second, the 
size of the median filter (BACK_FILTERSIZE) applied 
to the background grid used for smoothing large artifacts 
in the image. 

In Table |6l we list the completeness and reliability as 
a function of these parameters. Note that in this table, 
in order to reflect the true completeness for the speci- 
fied detection threshold, the completeness values of the 
SExtractor catalog for 6a, 7a and 10cr thresholds have 
been normalized by the percentage of sources in the DR2 
catalog which are beyond these thresholds respectively. 
Thus, the completeness at any threshold is expected to be 
100% if all the DR2 sources beyond that threshold are de- 
tected by SExtractor. The general trend observed is that, 
with increasing detection threshold, the completeness of 
SExtractor is fairly steady, but its reliability increases. 
However, for 5a and 6a detection thresholds, the com- 
pleteness is rather low (^75% or lower). Choosing higher 
detection thresholds such as 7-10cr does not seem to im- 
prove the completeness substantially. This is due several 
factors through which the source-finding in SExtractor 
seems to be different from the conventionally-used algo- 
rithms in radio astronomy. Firstly, for any given source, 
the reported peak flux density decreases as the detec- 
tion threshold is raised, although there is not much of a 
change in the reported rms. Hence, the SNR of sources 
close to the detection threshold decreases so as to be 
rejected by SExtractor. Choosing a 3a threshold usu- 
ally gives the correct peak flux density. This is also the 
reason why searching down to 3a and selecting sources 
at the desired, higher threshold, increases the complete- 
ness (although with the side effect of reduced reliability). 
The explanation for such an effect is hinted by the fact 
that the dimensions of the fitted sources decrease (this 
effect is quite significant for sources with SNR<10) with 
increasing detection threshold, implying that the num- 
ber of source pixels considered in the fitting process de- 
pends on the threshold. Secondly, several DR2 sources 
having SNR<10 and detected by other source-finding al- 
gorithms (excepting IMSAD), are not detected at all by 
SExtractor, not even at a reduced SNR. This might be 
due the differences in the fitting process rather than dis- 
crepant rms. Lastly, for some extended sources, posi- 
tions of the components reported by SExtractor are at 



least a few synthesized beams away from the positions 
given in the DR2 catalog; the latter almost always agree 
with the positions reported by other sources-findi ng al- 
gorithms (again, with the exception of IMSAD; see £ 4.1.4 
for anomalies associated with IMSAD). The completeness 
peaks at a background mesh-size of 5 beamwidths used 
with 3x smoothing, whereas reliability appears to in- 
crease steadily or remain constant with increasing mesh- 
size. BACKJF1LTERSIZE is not seen to change either 
the completeness or reliability much. SExtractor does 
not seem to achieve a high completeness and reliability 
(>85%) simultaneously for any given set of values for 
the input parameters. For reliability of > 90%, we see 
that completeness < 80%. Beyond a detection threshold 
of 6a, the reliability is quite good, yet less than 100% 
owing primarily to the acceptance of sidelobes and the 
fitting of extended sources handled differently. Never- 
theless, the highest reliability for SExtractor is better 
than the catalogs of other algorithms (except, maybe 
SAD), and hence it is best for cases where reliability is 
strongly favored over completeness. In such a case, we 
recommend setting the detection threshold to about 10ct, 
along with a large enough value for BACK.SIZE. T ests 
on SExtr actor with the ASK A P and IHancock et al.l sim- 
ulations (Huynh et al. 2011) suggest that the reliabil- 
ity is almost constant tor mesh-size between 10 and 100 
beamwidths, but increases with the detection threshold. 
They find that the completeness generally decreases or 
remains constant as mesh-size increases, but its change 
with respect to detection threshold is not explored. Our 
results are thus broadly consistent with the tests on the 
simulated images. 

The upper-left panel of Figure M shows the complete- 
ness of SExtractor sources as a function of their SNR 
in the DR2 catalog for mesh-size of 20 beamwidths and 
different values of the detection threshold. The com- 
pleteness for optimum input parameter values is shown 
in the upper panel of Figure [9] SExtractor approaches 
100% completeness faster than IMSAD, but slower than 
the other algorithms. This is possibly related to the ap- 
proach used for determining source pixels. The decrease 
in completeness beyond SNR~20 can be attributed to 
the extended-source components being reported differ- 
ently than other algorithms, as mentioned above. The 
lower panel of Figure [91 which shows the SExtractor re- 
liability in relation with other source-finding algorithms, 
reveals that this algorithm gives quite unreliable results 
for SNR<10. These low-SNR false sources are all point 
sources located near the edges of the DR2 image, where 
the rms is somewhat large (12 /iJy), but the rms reported 
by SExtractor is quite small (7-8 /iJy). This scenario is 
presented in the image cutouts in Figure 10 which also 
illustrate the ability of SExtractor to find sources with 
different morphologies located in different environments. 
Figure [TT] shows the SNR of sources detected by SEx- 
tractor, in comparison with the corresponding sources in 
the DR2 catalog. For SNR^30, the peak flux density 
of SExtractor sources sfind, SExtractor is generally less 
than that of the DR2 catalog sources, suggesting that 
the difference in the calculated peak flux density (and to 
a smaller extent, the associated uncertainty) is respon- 
sible for the observed departure of the SExtractor SNR 
from the DR2 SNR. 
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Fig. 8. — Completeness of the catalogs generated by different source-finding algorithms. Results for a background mesh-size (rmsbox) 
of 20 beamwidths , wh erever specifiable, are shown. For IMSAD and Aegean, the results are for the histogram option and csigma=l<T cmap 
respectively. See §4.1| for the definition of completeness used here. The decreasing completeness beyond SNR ~ 50 is due to deblending of 
multiple components of extended sources. The inset shows only the region where SNR is between 4 and 10. Smoothing^ over every 75 data 
points has been done before plotting. Note the low-number statistics for sources with SNR >70 as implied by Figure [ 



4.1.3. SAD (AIPS) 



algorithm finds all 
d (typically a mul- 



The "Search and Destroy" (SACF 
the pixels above a specified thresho 
tiple of the rms noise, which is assumed to be Gaus- 
sian) in the image, and merge contiguous pixels above 
the threshold into islands. The strength and size of each 
island is then estimated, followed by least squares Gaus- 
sian fitting of each island (if rms residual is too high, 
then multiple Gaussian fits may be applied). However, 
note that Gaussian statistics may not be a good model 
for the distribution of values in pixels well above zero 
flux density on accoun t of thermal noise, and calibration 
and imaging artifacts (ICotton & Peters 20111). The run- 
time for SAD searching tor sources down to bo in the DR2 
image is about 2 minutes. Additionally, to prepare the 
background rms image, the task RMSD takes ~15 minutes. 
For all the tests carried out on SAD, 31DEC11 AIPS was 
used. 



www.aips.nrao.edu/cgi-bin/ZXHLP2.PL?SAD 



SAD has several input parameters which affect the num- 
ber of sources detected. Here, we test the effect of the 
detection SNR (CP ARM), rms threshold, flux residual 
threshold (DP ARM, GAIN), and the size of rms box 
(IMSIZE). The parameter DP ARM (3) (along with 
GAIN added in quadrature) specifies the upper limit for 
the rms in the fitting box, while DPARM(7) and GAIN 
specify the upper limit on the residual flux in the fitting 
box. GAIN thus defines the fraction of the source flux 
which is acceptable in the residual image. By default, 
SAD uses the entire image to find the rms. We used the 
AIPS task RMSD to prepare rms images using mesh-sizes 
(specified by the IMSIZE parameter) of 5, 10, 20, and 
50 beamwidths. Although a decremental search in SNR 
(via CPARM) is recommended in the SAD help file, we 
found that such a search results in multiple sources being 
fit to a single genuine source during each iteration, espe- 
cially when the source is extended. Hence, we rejected 
this recommendation. As in the case of SExtractor, we 
searched for sources down to Ncr. Note that in SAD, the 
errors in the flux density are determined by theory from 
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the image rms (actnoise keyword in the image header). 

In Table [71 we list the completeness and reliability of 
the SAD catalog for different values of the input param- 
eters. We held the input parameters DPARM(3) and 
DP ARM (7) fixed at 1000 and 1 (in units of Jy beam" 1 ) 
respectively in order to get optimum completeness (with- 
out significant loss of reliability, as we found later, but 
possibly at the expense of correct SNR of the detected 
sources). Due to these large input values, we did not 
find any change in the results with the GAIN parame- 
ter. However, we found that reliability can be traded for 
completeness by setting smaller values for DP ARM (7) 
and GAIN. In general, the completeness and reliability 
of SAD increase with the detection threshold. Complete- 
ness increases with mcsh-sizc. So does reliability, though 
this quantity decreases significantly as the mesh-size is 
increased from 20 to 50 beamwidths, similar to sfind. 
Our inspection of the rms image for a mesh-size of 50 
beamwidths reveals that there are a few pockets where 
the rms is rather low (~1 /iJy). A profusion of false 
sources (or very low-SNR sources reported to have an 
SNR above the threshold) detected in these pockets are 
responsible for the markedly reduced reliability for the 
case of mesh-size equal to 50 beamwidths. Search with 
actnoise usually performs at least a few percent worse 
in terms of completeness and reliability than using mesh- 
size of 20 beamwidths. Curiously, the completeness and 
reliability of SAD is >90% for a wide range of input pa- 
rameters tested. We determine the optimum values of 
input parameters to be 10-20 beamwidths mesh-size and 
6-7ct detection threshold. Searching down to 3a followed 
by selection of sources greater than 6-7g may im prove 
completeness t o some extent. Huynh et al. (20111 and 
Hancock et al. (2012) have not tested SAD on simulated 



images. 

The completeness of SAD with respect to detection SNR 
and for different detection thresholds is shown in the 
upper-right panel of Figure[8j which depicts the high level 
of completeness close to the threshold and rapid increase 
with SNR compared to the other algorithms. The upper 
panel of Figure[9]plots the SAD completeness for optimum 
input parameters. Due to missed sources throughout the 
SNR spanned, the completeness is seen to hover close 
to unity, but not quite getting to 100%. The source- 
rejection criteria based on 2-D Gaussian fitting in SAD, 
defined by the several elements of the DP ARM input 
array, is the likely cause for such such missed sources. 
The lower panel of Figure [9] shows the reliability for the 
optimum values input parameters. As with the case of 
completeness, the reliability seems to hover close to unity. 
SAD provides the best reliability for sources SNR<20, be- 
yond which, sfind gives better results. Departure of the 
reliability from unity for a wide range of SNR values 
is due to the fitting of different components to extended 
sources. Some examples of how SAD fits different kinds of 
sources in the DR2 image, in c omparis on w ith other algo- 
rithms, is shown in Figure 10 Figure 11 shows the SNR 
of sources detected by SAD, in comparison with that of 
their counterparts in the DR2 catalog. Usually the SNRs 
agree with each other, but any disagreement between the 
two is due to the reported uncertainty in the peak. The 
peak flux densities reported by SAD match those in the 
DR2 catalog quite well. 

There is a possibility that the superior completeness 



and reliability of SAD may be due simply because the 
DR2 catalog, against which we are doing all of our com- 
parisons, was constructed from SAD. However, the con- 
struction of the DR2 catalog involved much more than 
running SAD (see ^3] for more details). In any case, an 
independent check using simulated sources is warranted. 
In Figure [l2j we compare the completeness and relia- 
bility of SAD w i th oth er source-finding algorithms on the 
Hancock et al. (2012) simulated image. We find a similar 
superior behavior of SAD. 




10 10' 10" 

Algorithm SNR 

Fig. 9. — The completeness (top) and reliability (bottom) of the 
catalogs generated by different source-finding algorithms. See j ]4.1| 
for the definitions of completeness and reliability used in this work. 
Results for a background mesh-size (rmsbox) of 20 beamwidths are 
shown. For IMSAD, these results are for the histogram option. For 
sfind, we have used «=10, whereas for SExtractor, SAD, IMSAD and 
Aegean, detection threshold of 3<r. Only those sources which are 
have then been selected from the respective catalogs prior 
to comparison with the DR2 catalog. The inset shows only the 
region where SNR is between 4 and 10. Smoothing over every 75 
data points has been done before plotting. Note the low-number 
statistics for sources with SNR >70 as implied by Figure [5] 



4.1.4. IMSAD (MIRIAD) 

IMSAEp^] computes the image rms noise by fitting a 
Gaussian to the image histogram, then searches for con- 

15 www.atnf.csiro.au/computing/software/miriad/doc/imsad.html 
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Fig. 10. — Examples of 5<r and brighter sources detected by various source-finding algorithms in the DR2 image. Top panel: 
blended components (30" cutouts centered on 03 h 32 m 32.2 s ,-28°03'09.4"), middle panel: source with sidelobes (20" cutouts centered 
on 03 h 32 m 06.1 s ,— 27°32'35.8"), bottom panel: region with a relatively large rms (12 /^Jy) at the corner of the image (3' cutouts centered 
on 03 h 31 m 19.4 s ,— 27°32'55.6"). The logarithmic flux density scale shown in each panel has units of Jy. The ellipses have major and minor 
axes and position angles according to the parameters reported by the respective algorithms. 
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Fig. 11. — The SNR of sources found by different source-finding 
algorithms relative to their SNR in the DR2 catalog. Smoothing 
over every 75 data points has been done before plotting. Note 
the low-number statistics for sources with SNR >70 as implied by 
Figure [2] 



tiguous pixels (islands) which are above some cutoff and 
fits the islands with Gaussian components. For fit- 
ting, the routine from the MIRIAD task imf it is used, 
whereas the island detection is based on the AIPS task 
SAD. IMSAD can also be used in a mode where the thresh- 
old background noise is user-specified. The run-time for 
this source-finding algorithm searching for sources down 
to 5er in the DR2 image is about 2 seconds. We used 
IMSAD version 8-May-98 from MIRIAD 4.2.3 (optimized 
for CARMA) to carry out these tests. 

As before, we searched for sources down to Ncr (spec- 
ified via the clip parameter), and then selected ^5er 
sources. The completeness and reliability were tested us- 
ing the histogram mode {noplt, hist options set) as well 
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Fig. 12. — Results of the completen ess (middle) and reliability 
(bottom) tests on the |Hancock et al.| simulated image. Matching 
radius of 15" was used to find counterparts. Input parameters to 
the source-finding algorithms are same as those given in Figure [9] 
Here, the co mple teness and reliability within each SNR bin are 
plotted (see j| 4.1t, which makes the se diagrams different from the 
ones given in Hancock et al.N2012t. For reference, the histogram 
of the simulated sources is also shown (top). The completeness at 
SNR~70 dropping below 90% is a result of highly-blended sources. 

as the user-specified mode (only the noplt option set; rms 
of 8 /xJy used) for background-rms determination. 

The results of these tests on IMSAD are given in Ta- 
ble [8] The noplt, hist mode usually performs better than 
the noplt mode at least by a few percent in terms of 
completeness, whereas the latter mode gives better com- 
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pleteness, again by a few percent. For achieving com- 
pleteness and reliability of 80% or more, the detection 
threshold needs to be ^6<r. Hence, the optimum use of 
IMSAD would be in this SNR cutoff regime, and the mode 
in which it is used will depend upon whether complete- 
ness or reliability is preferred. Searching down to 3c 
followed by rejection of sources below the desired thresh- 
old might give better completeness. The reason for the 
low completeness at a detection threshold of 5a is that 
several DR2 sources (which are detected by sf ind, SAD 
and Aegean) close to this threshold are not detected by 
IMSAD, not even at a reduced SNR. The reliability at 
this threshold is also rather low. This can be attributed 
to several false point sources detected near the edges of 
the DR2 image, which are mainly low-SNR features mis- 
taken to be high-SNR due to incorrect assumption of the 
local rms. Since IMSAD does not give the rms for each 
detected source, one needs to assume this quantity (a 
constant equal to 8 /xJy in our case), and thus, the SNR 
calculated in regions of increased rms, e.g. regions close 
to the edge of the DR2 image, is expected to be erro- 
neous. By using an rms image to find the local rms, this 
false-detection problem can be tackled. Another issue re- 
lated to false sources is that IMSAD does not reject a sin- 
gle sidelobe. All the sidelobes are reported as (genuine) 
sources. Thus, with respect to false sources detected in 
noisy regions, missed sources, and sidelobe recognition, 
IMSAD behaves like SExtractor. The sources detected by 
IMSAD in some of these scenarios are shown in the image 
cutouts in Figure 10 Hancock et al.| ( |2012[ ) tested IMSAD 
on a simulated image and found that its completeness 
and reliability is lower than other source-finding algo- 
rithms. However, through our tests on the DR2 image, 
it appears that IMSAD has a competitive performance if 
used with detection thresholds higher than ha. 

The completeness of the IMSAD catalog as a function 
of the detection SNR is shown in the lower left panel 
of Figure M The upper panel of Figure [9] shows the 
IMSAD completeness for optimum input parameters. Sev- 
eral sources having 10<SNR<20 as well as SNR>100, 
which are detected by other source- finding algorithms, 
are not detected by IMSAD. This is due to the inability 
of this algorithm to individually fit blended sources and 
components of extended sources; IMSAD tends to fit one 
elonga ted source for all components (see upper panel of 
Figure [To| . In Figure |9j the reliability of this algorithm is 
plotted lor the optimum input parameters. IMSAD shows 
reduced reliability between SNRs of 10 and 20, beyond 
which the reliability roughly flattens off at the 95% level. 
The reason for this reduced reliability is two- fold. First, 
several point sources detected close to the edge of the 
DR2 image have integrated flux densities less than their 
peak by a factor of a few, clearly indicating unphysical 
fitting of sources. These false sources can easily be re- 
jected by using a peak-to-integrated flux density ratio 
criterion. Second, IMSAD does not break up islands into 
components as mentioned above, resulting in extended 
source positions to be substantially different from those 
listed in the DR2 catalog. This was also found in the 



tests carried out by |Hancock et al. (|2012 1, which resulted 
in several false detections. Figure fTTl shows the SNR of 



of 

IMSAD sources compared with their counterparts in the 
DR2 catalog. The assumed constant rms as well as the 
somewhat lower peak flux density reported by IMSAD (for 



extended sources) are responsible for the lower SNR with 
respect to DR2 sources. 

4.1.5. Aegean 

^4 egearp'1 uses the FloodFill algorithm, which separates 
the foreground pixels from the background and groups 
them into islands. These "islands" are then passed on to 
the source characterization stage. Each island of pixels 
is fit with multiple Gaussian components. The number 
of components to be fit is determined from a surface- 
curvature map (cmap), derived from the input image 
with a Laplacian transform. Aegean thus performs a well- 
constrained multiple Gaussian fitting. A detailed de- 
sc ription of Aegean and its implementation can be found 
in jHancock et al. (2012), where the authors demonstrate 
this source-finding algorithm on a simulated image to 
produce catalogs for better reliability and completeness 
than other source finding algorithms. The run-time for 
Aegean searching for sources down to ha in the DR2 im- 
age is about 4 minutes (using two CPU cores). We used 
Aegean r706 (25-Jul-2012 release) for testing this source- 
finding algorithm. 

We tested the effect of the parameters: (i) innerclip 
(seedclip) and outer clip (floodclip), and (ii) csigma on 
the completeness and reliability of the Aegean catalog, 
innerclip defines the sigma clipping (lower limit) which 
is used for the detection ("seeding") of islands, whereas 
outerclip defines that used for reporting ( "flooding" ) the 
islands in the catalog, csigma is the sigma clipping 
parameter for the curvature map. Note that the r706 
version of Aegean uses an immutable mesh-size of 20 
beamwidths. We used innerclip — outerclip for our 
tests. Setting the innerclip to la lower than outerclip 
gives results alike equating these two clipping parame- 
ters, while using a la lower outerclip than innerclip seems 
to give a few percent better completeness and a slightly 
reduced reliability. 

Our results for Aegean are given in Table [9j csigma 
does not have a significant effect on either the complete- 
ness or reliability. As with our tests on other source- 
finding algorithms, we searched for sources down to Nc. 
The best completeness and reliability require a la or 
higher detection threshold. Using a 3cr detection thresh- 
old followed by the selection of greater than la might 
give better completeness at the cost of reliability Sur- 
prisingly, the completeness and reliability are not >95%, 
as seen in the tests with the Hancock et al. (2012) simu- 
lated image. This is because Aegean performs well with 
isolated point sources, but not quite with extended or 
blended ones. For extended sources, Aegean tends to 
fit some of the components at unexpected locations and 
with quite elongated Gaussians (the fitting-error flag is 
also set for these components) . Numerous small (but un- 
physical; size equal to the synthesized beam) components 
are also report e d for several extended sources. Since 
Hancock et al. (2012) have demonstrated that Aegean 
successfully fits all the components of blended simulated 
sources, we think that diffuse components and imag- 
ing artifacts such as negative sidelobes adversely affect 
the curvature map produced by this source-finding al- 
gorithm, which in turn determines the sources reported. 
An example of such a case is given in the upper panel 

16 www.physics.usyd.edu.au/%7Ehancock/index.php/Programs/Aegean 
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of Figure [10] Note that in the latest release of Aegean 
(r808; 13-Feb-2013) this situation is somewhat improved, 
since some of the elongated components from r706 are 
reported with more realistic source sizes (and also with- 
out any fitting-error flags) , and relatively fewer small-size 
components are fit. 



TABLE 5 

Completeness and Reliability of sf ind catalog 



rmsbox 




a 






(beamwidths) 


10 


5 2 


1 


0.1 


Completeness 


5 


80.9% 


77.8% 74.3% 


72.4% 


61.9% 


10 


85.7% 


83.7% 79.1% 


75.9% 


65.5% 


20 


88.4% 


84.5% 81.2% 


77.8% 


67.8% 


50 


89.1% 


85.7% 82.1% 


78.4% 


68.0% 


Reliability 


5 


78.8% 


84.9% 90.6% 


93.1% 


96.4% 


10 


86.8% 


91.6% 95.2% 


95.6% 


97.4% 


20 


88.7% 


92.5% 95.3% 


96.0% 


96.4% 


50 


87.5% 


91.4% 93.8% 


94.8% 


95.3% 



Notes— Only ^5<r sources have been selected from the 
SFIND catalogs prior to comparison with the DR2 catalog. 



TABLE G 

Completeness and Reliability of SExtractor catalog 



BACK.SIZE 




DETECT.THRESH 




(beamwidths) 


3(>5) 


5 


6 


7 


10 


Completeness* 


5 




82.2% 


72.1% 


76.8% 


79.1% 


77.9% 


5 + 


3 X smoothing 


83.4% 


73.7% 


78.7% 


80.0% 


79.0% 


5 + 


10 X smoothing 


83.3% 


73.0% 


78.4% 


80.0% 


79.0% 


10 




81.7% 


71.4% 


76.9% 


78.7% 


78.2% 


10 - 


1- 3 X smoothing 


82.1% 


71.1% 


76.7% 


79.1% 


77.6% 


10 - 


-f 10 X smoothing 


81.7% 


71.2% 


76.4% 


78.6% 


77.6% 


20 




81.6% 


70.6% 


75.9% 


78.6% 


77.6% 


50 




81.3% 


70.2% 


75.3% 


78.0% 


77.3% 


Reliability 


5 




70.4% 


87.5% 


95.1% 


96.3% 


97.0% 


5 + 


3 X smoothing 


72.4% 


89.0% 


96.1% 


97.2% 


97.4% 


5 + 


10 X smoothing 


72.2% 


89.5% 


96.1% 


97.2% 


97.4% 


10 




75.2% 


90.8% 


96.2% 


96.5% 


97.4% 


10 - 


+• 3 X smoothing 


75.8% 


90.1% 


96.2% 


97.0% 


97.4% 


10 - 


4- 10 X smoothing 


75.6% 


90.2% 


96.2% 


96.9% 


97.4% 


20 




77.0% 


90.3% 


96.4% 


96.9% 


97.4% 


50 




78.1% 


91.1% 


96.5% 


96.9% 


97.4% 



Figure [8] shows the completeness of the Aegean catalog 
as a function of SNR of sources in the DR2 catalog, and 
the upper panel of Figure [9] shows the completeness for 
csigma=la cmap , and searching for sources down to 3a. 
Aegean completeness appears to be quite good except for 
sources below 7a in the DR2 catalog, which are detected 
at a decreased SNR (below 5) and hence rejected from the 
catalog, and for a few components of extended sources 
above ~100<7, which are reported to have positions not 
matching those in the DR2 catalog. The decreased SNR 
reported for 5~7a sources is due to a slightly decreased 
peak flux density and a slightly increased rms with re- 
spect to the DR2 catalog. Figure [9] plots the reliability 
using csigma=la cmap , and searching for sources down to 
3a. The significant deviation of reliability from unity for 
SNR>10 sources, results from the several (false) small- 



TABLE 7 

Completeness and Reliability of SAD catalog 



IMSIZE 






CPARM 






(beamwidths) 


3(>5) 


5 


6 


7 


10 


Completeness* 


5 


77.8% 


75.5% 


80.0% 


83.7% 


78.4% 


10 


89.9% 


88.7% 


93.1% 


92.1% 


93.4% 


20 


93.4% 


92.2% 


95.9% 


95.4% 


97.1% 


50 


95.9% 


94.8% 


96.5% 


95.8% 


97.1% 


actnoise| 


87.4% 


86.7% 


95.2% 


96.7% 


95.7% 


Reliability 


5 


88.1% 


89.0% 


97.3% 


99.1% 


99.3% 


10 


91.2% 


93.0% 


97.4% 


97.4% 


96.6% 


20 


92.5% 


93.9% 


96.7% 


96.9% 


97.0% 


50 


85.1% 


87.3% 


88.1% 


88.6% 


92.1% 


actnoisej 


58.2% 


60.2% 


86.6% 


94.6% 


96.4% 



Notes— ^Search using the actnoise keyword in the FITS 
header (=7.465E-06 JY/BM). 



TABLE 8 

Completeness and Reliability of IMSAD catalog 



options 




clip 








3(>5) 


5 6 


7 


10 


Completeness* 


noplti. 


75.7% 


73.9% 78.1% 


84.2% 


79.3% 


noplt,hist 


75.4% 


76.5% 91.1% 


90.5% 


90.6% 


Reliability 


noplti. 


77.4% 


83.7% 95.3% 


96.5% 


96.8% 


noplt,hist 


78.2% 


79.2% 91.7% 


95.8% 


96.2% 



Notes— ^Clipping level is manually entered as the 
appropriate multiple of the background rms chosen 
to be 8 /iJy. 



TABLE 9 

Completeness and Reliability of Aegean catalog 



csigma 




inncrclip = 


outerclip 




(cmap rms) 


3(>5) 


5 6 


7 


10 


Completeness* 


0.5 


86.1% 


78.7% 83.8% 


86.1% 


85.8% 


1.0 


86.3% 


78.5% 84.2% 


86.4% 


85.8% 


2.0 


88.2% 


78.5% 84.5% 


86.8% 


85.5% 


Reliability 


0.5 


86.2% 


91.6% 91.6% 


92.1% 


91.7% 


1.0 


84.9% 


90.6% 90.5% 


90.5% 


91.1% 


2.0 


83.4% 


91.0% 91.2% 


91.4% 


89.9% 



Notes— *The fraction of sources in the DR2 catalog which 
arc ^6cr, 7cr and 10cr are 73.0%, 58.1% and 37.7% re- 
spectively. Completeness for these detection thresholds 
has been normalized accordingly. A detection threshold of 
3(>5) implies a search down to 3<r followed by the selection 
of only those source that are greater than 5o\ 



size, and elongated components rep orte d for extended 
sources, as mentioned above. Figure [TT] shows the SNR 
of sources detected by Aegean, in comparison with the 
corresponding sources in the DR2 catalog. The disagree- 
ment in SNR for h-la sources in the DR2 catalog has 
been explained above. For SNR>200 sources, the peak 
flux densities agree quite well between the Aegean and 
DR2 catalogs, but the rms reported by Aegean is con- 
sistently higher, which causes the disagreement at the 
high-SNR end of the diagram. 
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4.1.6. Summary of results from the efficiency tests 

We find remarkable differences between algorithms in 
terms of components fitted for extended sources, sidelobe 
rejection, and point sources detected in regions where the 
rms is appreciably larger than the mean rms. Our re- 
sults for c ompleteness and rel iabil ity are broadly simila r 
to those of |Huynh et al. (2011 ) and Hancock et al. (2012 ). 
For applications that need both completeness as well as 
reliability, sf ind and Aegean are good. Additionally, we 
found that the SAD algorithm within the widely avail- 
able AIPS package had a better performance. IMSAD 
also gives a good completeness and reliability for detec- 
tion thresholds ^6a. For transient searches, reliability 
takes preference over completeness, since false positives 
are likely to consume follow-up resources. Most tran- 
sient projects are likely to be searching in near real-time. 
However, in this particular project we were fortunate to 
have a deep reference image that was more than three 
times deeper than the single-epoch images. This allowed 
us to study reliability with real (rather than synthetic) 
datasets. From Figure |9j we see that the reliability 
of sfind is better than that of SExtractor, IMSAD and 
Aegean except for SNR near 8. However, the best relia- 
bility is provided by SAD. 



4.2. Transient Candidate Search 

Usin g the best-performing source finding algorithms 
from §4.1.6| and their optimum parameter values, we 
carried out a search for transient radio sources over all 
epochs. We ran sfind, SAD, IMSAD and Aegean on single- 
epoch images and obtained 49 single-epoch catalogs for 
each algorithm. We required that any potential transient 
candidate identified in the single-epoch catalog obey the 
following constraints. 

1. The source is not found in the reference catalog 
(within 2") of persistent sources. The reference 
catalog was constructed similar to DR2 but the 
sources were selected up to the 20% power point 
of the beam (i.e. 8 < 21.5' radius) of each pointing 
rather than the 34' interior region shown in Fig- 
ured] 

2. The source is at least a 7a detection. 

3. Is a genuine point-like source, i.e. it has 

(a) 0.9 < Sintl ! Speak < 1-5 

(b) a < 2.8"x 1.5, b < 1.6" X 1.5 

(a =major axis, b =minor axis; recall that the 
synthesized beam is 2.8" x 1.6") 

4. The source is at least 20 synthesized beams (20 x 
geometric mean of FWHMs; 42") away from the 
nearest: 

(a) bright source (> 500 jtiJy), so that any side- 
lobe emission is rejected 

(b) extended source 

The multiplicative factor of three-halves used in the 
selection of point sources, as well as the distance of 20 
beamwidths used for constraining the proximity from 
bright and extended sources, is somewhat arbitrary, but 
is based on several iterations of our transient search 



code and inspection of the cutouts of the resulting tran- 
sient candidates. By investigating how the major axes 
of sources increase with their distance from the point- 
ing center, we found that 9 ~ 21.5' appears also to be 
the threshold beyond which bandwidth smearing coupled 
with our constraints on the major (and minor) axes start 
rejecting genuine point sources. 

Due to the large number of synthesized beams searched 
(n = 1.8x 10 7 ) in this dataset, there is a modest probabil- 
ity that a transient candidate is due to nois 

< fH We thus 

carried out an analysis similar to Frail et al. (2012) (see 
Appendix A of that paper) to determine the SNR above 
which the probability of having the highest value of n 
Gaussian random numbers is ^1%. This corresponds to 
an SNR of 6.1 . How ever, following the recommendation 
oflFrail et all (120121) to have a higher SNR cutoff when 



the noise was not strictly Gaussian, we chose 7a as the 
lower limit for finding transients. 

The search method outlined above may miss transients 
which arc bright enough to be present in the reference 
catalog. Therefore, we also searched (with similar con- 
straints as above) for sources which are detected in the 
reference catalog of persistent sources and detected in 
only one of the single-epoch catalogs. 

For sfind, we adopted the parameters a=10 and 
rmsbox=20 for transient search. We found five candi- 
dates which are present only in a single epoch above the 
7a detection threshold. They also show up in the ref- 
erence catalog. However, all of these candidates seem 
to be variable but persistent, and are detected at a low 
SNR (~3-5c) in other epochs. Thus we do not find any 
transient with sfind. 

For Aegean, we used a clipping level of 7a and a 
curvature-map cutoff of la. We discarded all the tran- 
sients corresponding to islands which were too small to 
give a 6-parameter Gaussian fit (sources with flag 10000 
or 00100 set) since otherwise we were dealing with a large 
number of transients. Using these constraints, we found 
one candidate which is present in the reference catalog 
and is a persistent source detected at a largely reduced 
SNR in other epochs. Accordingly, Aegean does not yield 
any transients. 

For SAD, we searched for sources down to 3a with an 
rmsbox of 20 beamwidths, and used the input param- 
eters DPARM(3),DPARM(7)=(1000,1) to ensure opti- 
mum completeness and reliability based on our efficiency 
test. We do not find any transients with SAD. 

With IMSAD, we searched for sources down to 7a with 
the histogram option set. The several sources found to be 
transient candidates, are only variables at our detection 
limit, as evident through visual inspection of the images. 
Thus, there are no transients reported by IMSAD. 

5. DISCUSSION AND CONCLUSIONS 

In this paper we have explored the time-domain prop- 
erties of a 1.4 GHz survey made toward the Extended 
Chandra Deep Field South (E-CDFS) region. Six mo- 
saic pointings toward E-CDFS were taken in 49 separate 
epochs over a period of three months. Single-epoch im- 
ages allow us to explore the transient and variable radio 
sky at sub-milliJansky levels on timescales of days, weeks 

17 Prom theory we know that the statistics of beam values of 
interferometric maps should follow a Gaussian distribution. 
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and months. We will now use these data to assess the 



degree of variability ( |5.1[ ) and the transient rate (f 5.2 1 
of the radio sky, an d pre dict what will be seen by future 
wide-field surveys ( |5.3[ ). 

5.1. Comparison of variability with previous surveys 
We found in our study that only a small fraction 
(7/599 = 1.2±p 7%) of the point sources in the E-CDFS 
showed any significant variability on day-week-month 
timcscales. Evidently, the sub-milliJansky radio sky at 
1.4 GHz is not highly variable. The only previous sub- 
mjy study at 1.4 GHz was from a single deep pointing 
toward the Lockman Hole. With sampling timescales of 



19 days, Carilli et al. ( 2003 ) found less than 2% of sources 
above 0.1 mJy to be highly variable. 

Our findings at sub-mJy levels are consistent with the 
several previou s 1.4 GHz studies at hig her flux den- 



sity thresholds. Thyagarajan et al. (2011b analyzed the 



8444 deg 2 of the FIRST survey and found only 0.5% 
(1627/279407) of sources above 1 mJy va ried signifi- 
cantly on timescales of minutes to years. |Frail et al.| 
( fl994| imaged a 2° region toward a gamma-ray burst 
on timescales of 1-96 days and found that fewer than 
~1% of the source s above a flux density of 3.5 mJy were 
strongly variable, de Vries et al. (20041 imaged a 120.2 
deg 2 area of Sloan Digital Sky Survey Stripe 82 finding 
<1.4% (123/9086) of the radio sources to be strongly 
variable (i.e. > 4g) above flux dens ities of 2 mJy on a 7 
year timescale. lOfek & Frail (2011 ) do a two-epoch com- 
parison of FIRST and IN VSS point sources brighter than 
5 mJy and find that only 0.1% (43/4367) vary by more 
than 4cr over timescale s ranging from about 300 to 1700 
days. Bannister et al. ( 2011a|b ) analyzed 22 years and 
2775 deg z of MOST observations at 0.84 GHz and found 
only 0.17% (53/29730) strong variables above 1 4 mJy on 
tim escales of days to minutes to years. Finally, |Croft et 



al 



(20111 used the Allen Telescope Array to survey a 
690 deg^ area at 1.4 GHz. They compared their catalog 
to the NVSS, finding that 0.1% (6/4408) of the sources 
were highly var iable on a timescale of 15 years. 

In Figure [i~3| we plot the differential source counts for 
the persistent radio sky at 1.4 GHz, normalized in the 
usual way by the Euclidean rate ( |Huynh et al. 2005 ) . The 
steep evolution of the AGN with decreasing flux density 
is apparent, as is the flattening of the source counts near 
1 mJy. The fractional variability appears to be low, at a 
level of one percent, among the sources greater than 100 
(tiJy. 

Comparing the variability properties at 1.4 GHz of 
our sub-mJy population with those at higher flux densi- 
ties, we find no obvious change in the fraction of strong 
variables. This is despite the fact that radio-loud AGN 
(which dominate at milliJansky levels) diminish in im- 
portance as radio-quiet AGN and star forming galaxies 
begin to populate the radio sky at lower flux densities, 
the latter c onstituting about 50% of the sub - milliJansky 



radio sky (iSmolcic et al. 2008 



Padovani 2011). On 



the other hand, in a sample of brightness temperature 
limited radio sources, it would be expected that the vari- 
ability would increase with decreasing flux density. 

Our optical identifications for the seven radio variables 
show a mix of AGN and star-forming galaxies. Six of 
these are consistent with a nuclear source based on the 



carefully matched astrometry. There is no evidence from 
this sub-mJy sample that we are seeing a new source of 
variable radio emission such as supernovae, gamma-ray 
bursts, etc. that would be expected to be offset from the 
nucleus. 

Radio variability appears to be a strong function of 
frequency. The best study to compare to this one is 
the 5 GHz survey of Ofek et al. (2011), since it used 
a similar observing cadence and the identical statistical 
measures of variability. In this case it was found that 
^30% of point sources brighter than 1.5 mJy at 5 GHz 
were variable. 

Some caution is warranted in comparing fractional 
variability between different experiments. Differences in 
the cadence, integration time, duration, and angular res- 
olution will have a tendency of reducing the variability 
amplitude. For example, the low angular resolution of 
some surveys may have the effect of reducing the strong 
variability from compact radio sources embedded in dif- 
fuse emission. This current survey with a synthesized 
beam of 2.8" x 1.6" has the highest angular resolutio n of 
any previous variability survey (see Ofek et aL||2011 1 . 

Another factor to consider is that the degree ol vari- 
ability which can be detected in a given experiment, as 
measured by the modulation index m, is a function of 
the signal-to-noise ratio. A source in our present survey 
would have to have a mean flux density of 1 mJy in or- 
der to detect a 10% modulation, while at the 5c limit of 
the DR2 catalog (~40 ^Jy) a source would have to have 
m ^ 1.5 in order to be identified as a significant variable. 
No strong variables (i.e. m > 0.5) were identified in our 
survey, but only the about 90 sources in the DR2 catalog 
are bright enough (^126 /xJy) to have been identified as 
a strong variable. In either case the fraction of signifi- 
cant or strong variables is less than a few percent of the 
sample. 

The most robust conclusion that can be drawn is that 
the variable radio sky at 1.4 GHz appears to be relatively 
quiet with only a fraction of a percent of sources vary- 
ing substantially over a wide range of flux densities and 
timescales. 

5.2. Limits on transient areal density and rate 

We searched our multi-epoch data for transients but 
found none. The search was conducted on each image out 
to a radius of 21.5' from the pointing center. The single- 
epoch area out to that radius is 0.40 deg 2 , or a total area 
of 20 deg 2 for all 49 epochs. However, the sensitivity of 
the VLA antennas is not uniform across this area. The 
primary beam response is well-described by a Gaussian 
with a half- width to half-maximum of 15', falling to 20% 
response at our search radius of 21.5'. At the pointing- 
center the 7a flux density limit was approximately 210 
/ijy for each epoch. 

In order to calculate a limit on the a real density of any 
putative transient population we follow Ofek et al. (2011 1 
and parameterize the source number-count function as a 
power law of the form 



k(> S) = K (S/S y 



(4) 



where S is the peak flux density, k(> S) is the sky surface 
density of sources brighter than S, kq is the sky surface 
density of sources brighter than So, and 7 is the power 
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S (mjy) 

Fig. 13. — Normalized 1.4 GHz differential radio source counts for persistent sources from |de"Z otti ct al. (2010) an d the norma lized areal 
density of transients (or limits) as a function of the flux density for various surveys at this frequency. [Bannister et al.| ( |2011a[ ) survey at 
0.84 GHz is colored differently than the other survey s. Most of the sur veys are upper limits (wed ge symbols) and the sampled phase space 
is shown by the grey shaded area. Upper limits from |Frail et al.| l |1994| ) and |Bower fc Saul| ( |2011^ do not explore any new part of the phase 
space (non-grey area), and hence have been left out or this diagra m^ Uur upper-lim it is labele d as "E-CDFS" . Three surv eys have transient 
detections so far, 2<r error bars for which are shown according to |Gehrels| (|1986| . Note that |Thyagarajan et all| l |2011^ and |Bannister e"t| 
a], | ( |2011a| ) may have identified a few strong variables as transients (see * ]5.2| l , which would make th eir detections move downward on this 
plot. Trie black solid line is the model for AGNs and star-forming galaxies from ICondon I (Tl984). Lines of constant areal density are 
shown as blue dotted lines. The ho rizontal dashed lin es are estimates for the areal density tor known and expected classes of long-duration 
radio transie nts taken directly fro m |Frail et aL] | |2012| . The areal density for Swift J1644+57-like tidal disruption events has been modified 
according to B erger et al. | ( [201 2[ ) to reflect tneir true rate at 1.4 GHz. Upper limits from the ASKAP-VAST surveys are estimated to be 
an order of magnitude or more below the rate of orphan gamma-ray burst afterglows, and having a rms sensitivity ranging between 10 fijy 
and 0.5 mjy. 



law index of the source number-count function. We as- 
sume, for simplicity, a homogeneous source distribution 
in Euclidean Universe so 7 = 3/2. The one-si ded 2-a 
uppe r limit on the areal density is thr ee events (|Gehrels| 
[1986| . Therefore, using Equation C5 in |Ofek et al.| | |l201l| r 
we find that the 2a upper limit on areal density to a flux 
limit of 210 /iJy is 18.0 deg -2 per epoch. Given that we 
have 49 epochs the 2a upper limit on the areal density 
is k(> 0.21 mjy) < 0.37 deg -2 . We can further estimate 
an upper limit on the transient rate assuming a dura- 
tion tdur less than the shortest time between epochs of 
5R(> 0.21 mJy) < 268(i dU r/0.5 day)" 1 deg" 2 yr" 1 . 

Our upper limit on the areal density of transient 
sources at sub-milliJansky levels can be compared with 
the predi ctions based on previous surveys. The |Bower et| 
al. ( 2007[ ) survey is a useful benchmark since their areal 
density dominates all known classes of transients. Adopt- 
ing their measured two epoch rate of k(> 0.37 mjy) = 
1.5 deg -2 and assuming a Euclidean source distribution 
(i.e. 7 = 3/2) we predict n(> 0.21 mjy) = 3.5deg~ 2 at 
the flux density limit of our current survey. 

An alternative way to look at our results is to compare 
our null detection to the expected number of Bower et 



al. transients ex pected in our dataset . We use the pa- 
rameterization of Fender & Bell (2011 1 for the predicted 



Bower et al. transient rate as a function of flux density. 



log 



deg" 



J - In,, -5.13 (5) 



where k is the snapshot rate, and S v denotes the detec- 
tion threshold of the observations at the pointing center 
(i.e. 7a = 210 /iJy). Integrating both sides of equation 
[5 over the azimuthal angle and in 9 out to 21.5' we get 
about 0.42 transients per epoch, if the Bower et al. ( 2007 1 
transients are real. Since we have 49 epochs, we expect 
to have about 21 Bower et al. transients in our E-CDFS 
dataset. 

Our search on the E-CDFS field suggests that the areal 
density of radio transients is an order of magnitud e or 
more below the rate measured by Bower et al. ( 2007 ) (i.e. 
<0.37 deg~ 2 vs. 3.5 deg -2 ). Alternatively, we find a 2a 
upper limit of <3 transients while the predicted number 
is ^21 transients. Our w ork therefore appears to support 
that of Frail et al. ( 2012 ) , which found from a re-analysis 
of the Bower et al. data, that the transient rate was as 
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much as an order of magnitude smaller than previously 
reported. This conclusion would be more robust if the 
spectral index of the putative tra nsient population was 



better known. Ofek et al. (20101 was able to use other 



surveys to const rain the spectra inde x a > (where S„ 



) . Since the |Bower et al. ( 2007 ) rates were derived 



based on observations made mostly at 5 GHz, our only 
data provide strong constraints for ^ a ^ 1.1. A 
population of op tically thick a ~ 5/2 sources with a rate 



similar to that of Bower et al. ( 2007 ) would be undetected 
in our 1.4 GHz E-CDFS fields. 

Many of the same variability surveys discussed in §5.1| 
were also sensitive to transient s. S uperposed on the ra- 
dio source-count plot of Figure [13] are the results of sev- 
eral of these transients surveys. Light grey shaded ar- 
eas represent the transient phase space covered by each 
of the surveys, white space represents open phase space 
for future narrow-deep or wide-shallow surveys. With 
few exceptions most of these transient surveys result 
in upper limits. Also shown are the normalized areal 
density of several known and expec ted classes of long - 
duration radio transients, based on Frail et al. ( 2012 ). 
The nominal rates for the putative |Bower et al. ( 2007 1 
sources are about two orders of magnitude above the 
tidal disruption o r Sw J1644+57-like objects. |Thyagara^1 
jan et al. (2011) define a transient as having either a 
single detection in the analyzed epochs, or the high- 
est flux density 5 times greater th an the next highest 
one (d etection/upper-limit). T hus, Thyagar ajan et al.| 
(2011) and 2-epoch surveys like |Banmster et al.| ( |2011a[ ) 
may identify a strong variable as a true transient , which 
will move the source-count of detected transients lower 
in Figure [13] 

5.3. Future Radio Surveys 

There are several facilities built or under construction 
that will be capable of synoptic imaging at 1.4 GHz. All 
of these facilities have the exploration of the time domain 
as part of their core science programs. The Australian 
Square Kilometer Array Pathfinder (ASKAP) and the 
Apertif instrument on the Westerbork Synthesis radio 
Telescope (WSRT) will be using focal plane array tech- 
nology to instantaneously image an instantaneo us field of 
view (FoV) of 30 deg 2 and 8 deg 2 , respectively (|Murphy| 
IJ2 
Me 



et al. 12012 Oosterloo et al. 



2()10|. Soui I) Africa, is ! ^rilt i- 



mg MeerKAT. an array ot 64 13.5 -m diameter dishes 



with a FoV of 1 deg 2 at 1.4 GHz ( |Booth et al.| 2009[ ) 
Finally, there is the newly refurbished Karl G. Jansky 
Array (VLA ) has 27 25-m dishe s with a FoV of 0.25 deg 2 
at 1.4 GHz ( jPerley et al][20lT ). 
We list the capabilities of each of these telescopes in Ta- 
Survey speed (SS), normalized here to the VLA, 



blc 10 



is a useful figure of merit for inter-comparison of survey 
capabilities of long duration transients and is expressed 
as 



SS oc BW x n(A e /T sys ) 2 



(6) 



where BW is the bandwidth, £1 the FoV, A e is the total 
collecting areas times the aperture e fficiency e P , an d T sys 
is the antenna system temperature (Cordes 2008). The 



relative survey speeds are only approximate since some 
of the system parameters for MeerKAT and ASKAP 
have not been confirmed with on-the-sky testing. Like- 



wise, we have assumed that radio frequency interference 
(RFI) limits the VLA bandwidth to only 50% of its max- 
imum BW. Larger fractions are achievable in the more 
extended array configurations and with better RFI exci- 
sion of the data. We have calculated the ASKAP survey 
speed with both the 18 phased array feeds which are cur- 
rently funded and with the full 36 as originally s peci fied. 
Despite these uncertainties, it is clear from Table [l0| that 
to within factors of a few, these are all powerful wid<> field 
imaging facilities. 



TABLE 10 
Telescope Specifications 



Telescope 


BW 


a 


D 


N 


e e 




SS 




(MHz) 


(deg 2 ) 


(m) 






(K) 




VLA 


512 


0.25 


25 


27 


0.5 


26 


1.0 


ASKAP 


300 


30 


12 


18(36) 


0.8 


50 


1.1(4.6) 


Apertif 


300 


8 


25 


13 


0.75 


70 


1.3 


MeerKAT 


750 


1 


13.5 


64 


0.7 


30 


4.1 



Notes— Here BW is bandwidth in MHz, f! is the field of view 
in deg 2 , D is the antenna diameter in meters, N is the number 
of antennas in the array, e e is the aperture efficiency, T 3ya is 
the system temperature in Kelvin, and SS is the survey speed 
normalized to the VLA. 



To illustrate these survey capabilities and compare 
them to what we currently know about the transient and 
variable radio sky, we will use the example of a electro- 
magnetic counterpart search for gravitational waves. For 
a good overview of the topic of EM- GW searches, and 
the m ain is sues, we refer the read er to |Metzger &: Berger 
( |2012[ ) and |Nissanke et al. | ( pOTTj ) . 

Long-duration radio emission has been predicted to 
originate from the merger of a neutron star binary from 



sev eral s ources including the merger shock (Kyutoku et 



al. 2012), afterglow emission from the bea med outflow 



m the relativ istic and non-relativistic phases ( Metzger & 



Berger 



2012), and from quasi-isotropic, mildly relativis 



tic outflows ejected dur ing the merger (Nakar & Piran 
2011[ |Piran et al.| 2012). All of these mechanisms de- 



pend on the amount ot energy put into shocked material 
and the density of the ambient medium. Predicted flux 
densities and timescales therefore vary over a wide range. 
We take as an example the detection of a signal with a 
peak flux density of 100 /iJy. Such a signal might be 
expected to occur for an ambient medium with a den- 
sity of 0.1 cm" 3 on day timescales for a merger shock, 
or on year-long timescales for mildly relativistic ejecta. 
We note that if short duration gamma-ray bursts are the 
dominant population of neutron star mergers, then the 
predicted radio signal would be difficult to detect with 
the telescopes in Table 10 given the the canonical e nergy 



and ambient density inferred for th is population (Fong 



et al.||2012| | Metzger fc Berger||2012[ ). 



1'he median sky localization ot a gravitational-wave 
source will be 60 deg 2 with a three-element GW net- 



work, and 7 deg 2 with a five- element GW network (Nis 
sanke et al.|2011 1. With such large error boxes, the mam 



challenge for the identification of a EM-GW counterpart 
will be to distinguish it from the foreground of false pos- 
itives. Optical-only searches for EM-GW counterparts 
are expected to be overwhelmed by false positives at the 
require d depth of 22-23 mag and special strategies are re- 
quired (Nissanke et al. 2011 ). However, as we have shown 
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from this paper, the radio sky at 1.4 GHz is relatively 
quiet. Integrating the differential source c ounts in Fig- 
ure 13 using the fit from Huynh et al. ( 2005 ), we estimate 
that the number of persistent radio sources above 100 fiJy 
to be 910 deg~ 2 . We have estimated that the fraction of 
strong variables, on a wide range of timescales, is likely 
to be one percent or less, or 9 strong variables per square 
degree. Further, the radio variables t hat we have seen to 
date have all been nuclear sources (j |3.3| ). Such variable 
or transient sources could be easily rejected as EM-GW 
counterparts since significant offsets are predicted from 
the host galaxy based on binary neutron-star popula- 
tion synthesis model s and measurements from short-hard 
gamm a-ray bursts (Belczynski et al. |2006 Fong et al. 
20101. Lastly, the number ot unrelated transients is also 
expected to be much less. Our derived limit on the tran- 
sient areal rate of k(> 0.21 m Jy) < 0.37 deg -2 translates 
to a <1.1 deg~ 2 at 100 fiJy (f|5.2|) for a Euclidean distri- 
bution. 



The l imit on the know n transient populations is 
even smaller ( |Frail et al.||2012[ ). 

We believe that the multi-wavelength approach which 
we have taken here should inform future searches. For 



the radio variables that we found in the E-CDFS, we were 
able to identify the s ource of the emission using optical 
images and spectra (j |3.3[ ). A similar strategy could be 
employed to identify lalse-positives for the small num- 
ber of radio variables or transient sources identified in 
EM-GW counterpart searches. Whether this approach 
will ultimately lead to a robust EM-GW counterpart de- 
tection is uncertain, but in terms of characterizing the 
variables, minimizing false-positives and getting an early 
sense on the nature of transients, we suggest that joint 
radio-optical searches will be fruitful for exploring the 
dynamic sky. 
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